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ABSTRACT 

An  analytical  and  experimental  investigation  of  the  inertial  stresses 
resulting  from  very  high  internal  heat  generation  rates  is  described.  The 
results  show  that  the  amplitude  of  the  thermally  induced  elastic  vibrations 
in  an  unrestrained  thin  rod  are  dependent  upon  the  ratio  of  heating  time  to 
the  fundamental  period  of  oscillation,  and  that  the  most  severe  stress 
^amplitude  results  from  a  constant  heating  rate  as  compared  to  an  exponential 
type  of  input. 

In  the  experimental  investigation,  rapid  heating  was  accomplished  by 
the  discharge  of  a  low- inductance  capacitor  bank  (0.1  uh,  2800  Joules) 
through  the  rod.  A  modified  biref ringent-coating  method  was  used  for 
observation  of  the  thermally  induced  oscillations.  Utilizing  the  frozen 
stress  technique  of  three-dimensional  photoelasticity,  an  initial  pattern 
of  interference  fringes  was  introduced  into  thin  strips  of  birefringent 
material  which  were  securely  bonded  to  each  side  of  the  rod.  Strain 
oscillations  were  measured  by  direct  observation  of  the  movement  of  the 
fringe  pattern  with  a  high  intensity  incident  light  source,  a  crossed- 
circular  polariscope,  and  a  high  speed  framing  camera.  Interpretation 
of  these  measurements  required  a  dynamic  calibration  and  a  correction 
factor  for  the  reinforcing  effect  of  the  coating.  Oscillations  with 
periods  from  35  to  250  microseconds  and  amplitudes  of  9000  psi  were 
measured.  The  results  were  In  good  agreement  with  the  uncoupled  thermo- 
elastic  theory,  and  it  is  concluded  that  this  technique  is  suitable  for 
transient  measurements  in  the  presence  of  large  magnetic  fields  which 
normally  restrict  the  use  of  electronic  methods. 
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1.  INTRODUCTION 

The  effect  of  high  heating  rates  on  the  behavior  of  mechanical 
systems  is  becoming  a  problem  of  greater  and  greater  engineering 
importance.  For  example,  the  increased  use  of  nuclear  reactors  for 
power  generation  and  for  studies  of  intense  radiation  damage  to 
materials  has  resulted  in  a  need  for  an  understanding  of  the  effects 
of  rapid  internal  heat  generation  on  the  mechanical  response  of 
reactor  components.  In  a  fast  burst  reactor  the  sudden  criticality 
results  in  a  very  high  fission  rate  accompanied  by  a  rapid  increase 
in  temperature  of  the  reactive  components.  As  a  result  of  the 
inertia  of  these  components  rapid  thermal  expansion  cannot  take 
place  and  large  transient  inertial  stresses  are  produced.  Hence, 
the  designer  is  faced  not  only  with  the  problems  of  temperature 
gradients  upon  cooling  and  the  resulting  quasi-static  thermal  and 
contact  stresses,  but  also  with  initial  inertia  effects  which 
govern  the  structural  response  during  the  rapid  temperature 
change. 

In  order  for  significant  inertial  stresses  to  be  developed, 
it  has  been  found  that  the  imposed  heating  rates  must  cause  a 
substantial  temperature  change  in  a  structure  in  times  short 
compared  to  the  mechanical  response  time  of  the  system.  For 
large  thermal  fluxes  incident  on  a  surface  or  for  rapid  internal 
heat  generation  rates,  these  temperature  changes  can  occur  in 
times  which  are  short  relative  to  the  mechanical  response  time 
of  the  heated  structure.  Hence,  the  ratio  of  the  time  of  heat 
input  to  a  characteristic  mechanical  time  will  be  indicative  of 
the  severity  of  the  inertial  stresses  caused  by  rapid  heating. 
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The  mechanical  response  times,  or  periods  of  free  vibration,  for 

structural  members  such  as  beams,  plates,  rods  and  shells  usually  vary 
-4  -1 

anywhere  from  10  to  10  seconds;  consequently,  inertial  stress 

effects  become  important  for  heating  times  which  may  vary  from  less 
•6  -2 

than  10  seconds  to  as  high  as  10  seconds.  It  is  evident,  then, 
that  even  though  the  heating  time  may  be  long  in  an  absolute  sense, 
inertial  effects  may  still  create  serious  design  problems.  In  addition, 
the  problem  of  coupling  effects  (i.e.,  when  the  temperature  field  is 
significantly  affected  by  the  state  of  deformation  of  a  body)  between 
the  elastic  and  thermal  field  equations  may  exist  when  the  rate  of 
strain  becomes  high.  However,  both  the  inertia  and  coupling  effects 
may,  within  the  bounds  of  practicality,  be  treated  separately,  and 
for  the  present  study  it  will  be  shown  that  coupling  effects  are 
negligible  and  attention  can  be  restricted  to  inertial  effects. 

A  literature  review  indicates  that  only  within  about  the  last  ten  \ 
years  has  considerable  attention  been  given  to  these  problems.  Actually, 
as  early  as  1837,  in  the  very  beginning  of  the  development  of  thermo¬ 
elasticity,  Duhamel^  considered  inertial  effects  by  noting  that  if  the 
temperature  rate  is  slow  enough,  inertia  terms  in  the  field  equations 
of  thermoelasticity  can  safely  be  neglected.  However,  as  far  as  can 

be  determined,  the  first  extended  investigation  of  inertial  effects 

2 

was  carried  out  in  1950  by  Danilovskaya  for  the  case  of  an  elastic 

half-space  suddenly  heated  at  the  boundary.  Further  work  was  done 
3 

in  1952  by  Mura  who  determined  the  transient  stress  distributions 
in  an  elastic  half-space  and  an  infinitely  long  solid  cylinder  sub¬ 
jected  to  a  sudden  temperature  change  at  the  surface.  The  results 
show  that  in  addition  to  the  thermal  expansion  wave,  an  elastic  stress 
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wave  arising  from  the  inertial  effect  is  propagated  through  the  half¬ 
space  at  the  sonic  velocity,  while  for  the  cylinder,  radial  vibrations 
occur  but  are  of  extremely  small  amplitude  even  for  the  most  severe 
type  of  heating  -  a  step  change  in  surface  temperature.  Additional 

work  on  the  inertia  effects  for  the  elastic  half -space  was  later 
4  5 

performed  by  Nowacki  and  Ignaczak  in  1957  and  by  Sternberg  and 
£ 

Chakravorty  in  1958.  All  of  this  work,  however,  had  no  direct 
practical  application,  but  it  did  serve  as  a  valuable  aid  to  the 
understanding  of  the  inertia  effects  associated  with  rapid  heating 
of  an  elastic  body. 

Space  exploration  activities  directed  attention  to  the  serious 
heating  problems  encountered  in  the  design  of  high-speed  propulsion 
units,  particularly  aerodynamic  heating  of  reentry  vehicles.  In  1955 
Boley7  analyzed  the  behavior  of  a  simply  supported  rectangular  beam 
of  depth,  h,  and  thermal  diffusivity,  D,  subjected  to  a  sudden 
application  of  heat  to  the  upper  edge.  He  used  the  elementary 
Bernoulli-Euler  theory  of  bending  with  the  assumption  that  the 
temperature  remains  constant  along  the  beam  and  varies  only  through 
the  depth.  Both  these  assumptions  gave  good  practical  approximations 
to  the  aerodynamic  heating  of  the  leading  edge  of  wing  sections.  His 
calculations  show  that  the  inertial  effects  become  important  when  a 
characteristic  heating  time  defined  by  h  /d  is  less  than  half  the 
period  of  free  lateral  vibration.  As  the  heating  time  approaches  zero, 
the  maximum  dynamic  deflection  approaches  twice  the  static  deflection 
which  would  result  if  the  beam  were  heated  slowly  along  the  upper  edge, 
and  as  the  heating  time  approaches  infinity,  the  inertia  forces  dis¬ 
appear  and  the  static  solution  alone  remains.  A  year  later  Boley  and 
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Barber  extended  this  work  to  include  the  analysis  of  thermally  induced 
vibrations  of  rectangular  plates  and  beams  for  several  typical  heat 
applications  encountered  in  spacecraft  design. 

Recent  developments  in  reactor  technology  have  resulted  in  consider¬ 
able  thermal  stress  analysis  of  components  subjected  to  internal  heat 

generation.  The  most  significant  quasi-static  problems  investigated 

9 

include  the  work  of  Kammash,  Murch  and  Naghdi  (1959),  who  worked  out 

the  complete  elastic-plastic  solution  for  the  cylinder  subjected  to  a 

radially  distributed  heat  source,  and  the  work  of  Weil^  (1962),  who 

determined  the  thermal  stresses  in  cylinders  subjected  to  internal  heat 

generation  due  to  7-ray  radiation.  The  only  investigation  known  to  the 

author  of  inertia  effects  associated  with  rapid  internal  heat  generation 

has  been  performed  by  Burgreen*1  (1962), who  has  determined  the  dynamic 

elastic  response  of  fuel  elements  having  geometries  of  thin  rods,  thin 

spherical  shells,  and  solid  spheres.  His  calculations  reveal  the  nature 

of  the  inertial  stress  wave  propagation,  and  he  concludes  that  the 

maximum  amplitude  of  these  waves  is  governed  by  the  ratio  of  the  heating 

time  to  the  fundamental  period  of  free  vibration  of  the  element.  These 

7  8 

results  agree  in  principle  with  those  obtained  by  Boley.  ' 

In  all  of  the  above  work  mentioned  thus  far,  the  thermo-mechanical 

coupling  between  the  elastic  and  thermal  energy  equation  derived  origin- 
12 

ally  by  Jeffreys  there  is  an  additional  energy  term  which  accounts  for 
the  work  done  by  the  volume  change  and  manifests  itself  as  a  thermo- 
elastic  dissipation  which  is  an  important  cause  of  internal  damping. 
Because  of  the  heating  or  cooling  associated  with  a  change  of  volume, 
the  temperature  distribution  is  dependent  upon  the  state  of  deformation, 
and  hence,  the  exact  solution  to  a  dynamic  thermoelastic  problem  is 
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complicated  by  this  coupling  between  the  dynamical  equilibrium  equations 

and  the  general  energy  equation,  or  heat  conduction  equation. 

Recently,  much  attention  has  been  given  to  the  determination  of 

such  rigorous  solutions.  Fundamental  to  this  effort  has  been  the  work 
13 

of  Biot  who  has  presented  a  unified  treatment  of  thermoelasticity  by 

application  and  extended  development  of  the  methods  of  irreversible 

thermodynamics.  This  work  shows  that  for  time  dependent  problems,  the 

effects  of  the  thermal  state  and  the  deformation  of  an  elastic  solid 

cannot  logically  be  treated  separately.  Thus  far,  only  one  partial 

rigorous  solution  of  this  nature  has  been  published,  but  some  effort 

has  been  directed  toward  solution  of  the  coupled  equations  which  are 

harmonic  in  the  time  variable  and  one  spatial  coordinate.  In  1957 

investigations  of  plane  waves  in  a  thermoelastic  solid  were  carried 

out  by  Deresiewicz, Lessen,^5  and  Chadwick  and  Sneddon.1^  These 

authors  agree  that  even  for  metals  at  reasonably  high  temperatures 

the  effect  of  coupling  between  elastic  and  thermal  motions  can  be 

neglected.  Lockett^  (1957)  has  also  worked  out  the  two-dimensional 

solution  of  the  propagation  of  Rayleigh  waves  and  found  that  coupling 

effects  introduce  a  difference  of  less  than  one  percent  in  the  velocity 

18 

and  amplitude  of  these  waves.  Finally,  in  I960  Chadwick  has  provided 

an  excellent  survey  of  the  present  dynamical  theory  of  thermoelasticity 

in  which  he  discusses  the  irreversible  thermodynamics  of  an  elastic 

solid,  provides  a  derivation  of  the  general  thermoelastic  equations, 

and  examines  the  characteristics  of  the  general  thermoelastic  boundary 

value  problem.  He  applies  his  results  to  some  specific  problems  which 

nicely  illustrate  the  nature  of  thermoelastic  dynamics.  In  their  text, 
19 

Boley  and  Weiner  have  summarized  much  of  the  above  work  which  took 
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place  before  1959,  and  have  worked  out  numerical  results  to  specific 
problems  in  order  to  illustrate  the  effects  of  inertia,  coupling  and 
thermoelastic  damping.  They  have  demonstrated  that  for  most  engineering 
design  purposes  the  elastic  and  thermal  field  equations  may  be  uncoupled, 
which  will  allow  considerable  simplification  of  solutions  to  dynamic 
problems  in  thermoelasticity. 

All  of  the  aformentioned  work  has  been  analytical,  but  experimental 
verification  of  the  results  is  noticeably  absent.  Despite  the  consider¬ 
able  number  of  studies  which  have  indicated  that  rapid  heating  of 
structures  is  accompanied  by  serious  Inertial  stresses,  to  this  author's 

knowledge,  there  have  been  only  two  experimental  investigations  to  verify 

20 

these  results.  In  1961,  Michaels  investigated  the  effect  of  incident 
thermal  ratlation  on  one  face  of  a  cylindrical  rod  and  successfully 

21 

demonstrated  the  existence  of  the  resulting  elastic  wave.  Llndholm 
(1962)  designed  an  arc  image  furnace  which  was  used  to  produce  an 
extremely  high,  short-duration  thermal  flux  on  the  surface  of  a  thin 
spherical  cap.  His  results  consisted  primarily  of  subjective  observa¬ 
tions  of  the  gross  effects  of  the  rapid  heating,  and  were  inconclusive. 

In  view  of  this  situation,  it  was  concluded  that  additional  experimental 
work  was  highly  desirable.  Realizing  that  serious  practical  problems 
surround  the  design  and  operation  of  suitable  experimental  techniques, 
considerable  thought  was  given  to  the  approach  to  be  taken.  This  led 
to  the  decision  to  conduct  a  study  of  the  behavior  of  a  long,  thin  rod 
subjected  to  high  internal  heat  generation  rates.  Uniform  heating  was 
accomplished  electrically  by  the  discharge  of  a  low- inductance  capacitor 
bank  through  the  rod.  The  movement  of  the  isochromatic  fringe  pattern 
in  a  photoelastic  coating  on  the  bar  was  observed  with  a  high  speed 
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framing  camera.  The  significant  feature  of  this  optical  method  is  that 
the  transient  axial  strain  can  be  observed  during  the  heating  pulse. 
Semiconductor  strain  gages  on  the  bar  were  also  used  to  observe  the 
steady  state  oscillations  which  occurred  after  the  electrical  heating 
pulse  and  to  provide  a  dynamic  calibration  of  the  photoelastic  coating. 
In  preparation  for  these  experimental  techniques  a  one-dimensional 
analysis  of  the  elastic  behavior  of  the  bar  was  carried  out  to  provide 
a  theoretical  comparison  for  the  experimental  measurements.  The  results 
of  these  activities  are  shown  and  discussed  in  the  following  sections. 


* 


8. 


2.  THEORETICAL  ANALYSIS  OF  THE  TRANSIENT  ELASTIC 
STRESS  DISTRIBUTION  IN  A  ROD  SUBJECTED  TO 
RAPID  INTERNAL  HEAT  GENERATION 

The  system  selected  for  study  consisted  of  a  long,  completely 
unrestrained  elastic  bar  of  length  2£  and  constant  cross-sectional 
area  A  (see  Fig.  1)  subjected  to  a  uniform  internal  heat  generation 
rate  per  unit  volume,  q.  If  the  heating  is  rapid,  the  inertia  of  the 
bar  will  oppose  thermal  expansion  and  stresses  will  develop.  Also,  if 
the  heating  time  is  short  relative  to  the  time  required  for  heat  transfer 
away  from  the  bar,  the  temperature  will  be  uniform  and  independent  of 
position  along  the  bar.  It  will  be  assumed  that  only  simple  tension  or 
compression  exists  and  that  the  axial  stress  is  accompanied  only  by 
lateral  strains  -  vex»  This  assumption  is  valid  provided  the  wave¬ 
length  of  axial  waves  is  large  compared  to  the  lateral  dimensions  of 
22  23  24  25 

the  bar.  *  '  ’  Hence,  inertial  effects  corresponding  to  lateral 
motion  will  be  neglected  so  that  the  following  analysis  is  restricted 
to  a  one -dimensional  problem  which  further  implies  that  plane  sections 
remain  plane  during  deformation. 

Let  the  origin  of  a  set  of  rectangular  Cartesian  coordinates  (x,y,z) 

coincide  with  the  center  of  a  bar  such  that  the  x-axis  lies  along  the 

axis  of  the  bar.  For  one -dimensional  motion  of  the  bar,  only  the  axial 

26 

stress  o^  remains,  and  the  stress  equations  of  equilibrium 


do 

XX 

“sr  + 

do 

tF  * 

So 

xz 

“ST  " 

P  o 
dtz 

do 

do 

do 

-,2 

-aF-  * 

tF  + 

tF  ■ 

.  o  v 
P  2 
dtz 

do 

do 

do 

.2 

"sf  + 

tF  + 

zz 

“sr  " 

_  a  w 

P  — f 
dtZ 

(2.1) 


9. 


reduce  to 


(2.2) 


Since  this  work  will  consider  only  elastic  deformation,  the  behavior  of 
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the  bar  Is  governed  by  the  linear  thermoelastic  stress-strain  relations 


°x  "  X(€x  +  cy  +  €z)  +  2G  ex  -  (3k  +  2G)  a  (T  -  Tq) 

0y  =  X(€x  +  ey  +  ez)  +  2G  €y  -  (3X  +  2G)  a  (T  -  Tq)  -  0  (2.3) 

a  -  X(cv  +  e  +  e  )  +  2G  «  -  (3X  +  2g)  a  (T  -  T  )  -  0  . 

z  x  y  z  z  o 

where  the  temperature  Tq  corresponding  to  the  stress-free  state  has 

been  added,  and  all  shear  components  of  stress  and  strain  are  zero. 

Eliminating  e  and  6  ,  Eqs.  (2.3)  reduce  to 
x  y 

ax  "  E  tex  '  a(T  '  VJ  5  (2,4) 

also, 

«y  -  «8  -  -  vex  +  (1  +  v)  a  (T  -  t0)  .  (2.5) 
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From  the  strain-displacement  relations  for  small  strains, 


(2.6) 


€ 

y 


e 

z 


dw 

57  * 


(2.7) 


Since  the  surface  tractions  are  specified  ( a ■  0  at  x  ■  i£) ,  It 
will  be  convenient  to  write  the  equation  of  motion  (2.2)  in  terms  of  the 
stress.  Hence,  differentiating  (2.2)  with  respect  to  x  and  substituting 
(2.6)  into  (2.4)  and  then  combining  with  the  partial  derivative  of  (2.2) 


i 


10. 


we  gee  the  one-dimensional  Beltrami -Michell  equation  of  compatibility 


a2a 

_ _ 2 

Sx2 


£ 

E 


a2* 
_ _ > 

at2 


p  a 


a2T 

at2 


(2.8) 


If  the  temperature  history  of  the  bar  is  known  or  prescribed, 
solution  of  Eq.  (2.8)  will  yield  the  stress  in  the  bar  directly.  However, 
in  most  problems  the  heat  generation  function  is  usually  given  and  the 
temperature  calculated  from  the  energy  equation.  For  a  linear  isotropic 

solid,  the  complete  energy  equation  is  well  known.  It  was  obtained 
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originally  by  Jeffreys  and  later  re-derived  and  discussed  by  Biot. 

This  equation  is  sometimes  referred  to  as  the  coupled  heat,  or  diffusion 

equation,  and  is 


k  v2T  +  q  -  7cv  +  (3\  +  2G)  a  Tq 


(2.9) 


where  A  **  ■  elastic  dilatation.  The  addition  of  the  last 

term  in  the  heat  equation  accounts  for  the  work  done  by  deformation  of 
the  solid  at  temperature  Tq.  Since  it  has  been  assumed  that  the  temper¬ 
ature  is  uniform,  the  first  term  on  the  left-hand  side  is  identically 
zero.  From  the  thermoelastic  stress-strain  relations 


Sa 

ST 


3\  +  2G 


s®  . 

3I  + 


H 

(2.10) 


8  -  o  +  o  +  a 
x  y  z 

and  for  the  one  dimensional  formulation  8  =  o^.  Substituting  in  (2.9) 
and  differentiating  once  with  respect  to  time 
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5? 


■)  n2  b2a 

[7c  +  3a  T  (3\  +  2G)]  ~  +  a  T 
v  °  bt2  0  bt2 


(2.11) 


which  when  combined  with  (2.8)  yields 


32a 


dx 


x  p 
2  "  E 


_  2 
r  Ea  i 

,  o 

^2 

pa 

1  2 
7cv+3(3\+2G)a  Tq 

3t2  * 

Pcv+3(3X  +2G)a2TQ 

s?  .  <2-12> 


The  problem  has  obviously  been  considerably  simplified  by  the 
assumption  that  the  temperature  is  independent  of  x.  Otherwise  Eqs. 
(2.8)  and  (2.11)  would  be  coupled  by  the  last  term  in  (2.11),  and  would 
require  the  simultaneous  determination  of  temperature  and  stress.  This 
simplifying  assumption  is,  however,  justified  (see  Appendix,  Sect.  2). 
Equation  (2.12)  can  be  simplified  further  by  noting  that  for  most  metals 


JEa2T 


7c  +3 ( 3A+2G ) JO2 T 


<  <  1 


> 


and 

pcv  »  3(3\+2G)JaZTo  , 

where  J  *  1.07  x  10  4  BTU/in-lb  has  been  inserted  for  dimensional 
correction  where  appropriate.  For  instance,  for  304  stainless  steel 


7cv  -  (0.28) (0.12)  -  0.034  BTU/ln3cF 

2 

3(3X+2G)Ja2TQ  -  3(73xl06)(10'4)(9xl0"6)  (70)  -  0.00016  BTO/in3°F 
JEa2T 

- 2 - =—  -  0.0006  3 

7cv+3(3X+2G)Ja  Tq 


Hence,  the  above  inequalities  are  valid 


12. 


It  is  interesting  to  note  that  the  bracketed  coefficients  in  the 

second  and  last  terms  of  Eq.  (2.12)  can  be  interpreted  as  the  corrections 

for  the  adiabatic  sonic  velocity  and  the  specific  heat  at  constant 

pressure,  respectively.  These  expressions  agree  with  the  general  three- 
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dimensional  derivations  given  by  Jeffreys,  namely  that 


adiabatic 


E 

P 


1  - 


EC^TqJ 


7cv+3(3\+2G)Ja  Tq 


where 


isothermal 


JJ 


sonic  velocity. 


(2.13) 


Hence  for  this  work 


3ad.  *  aisoth. 


(2.14) 


Also,  the  specific  heat  at  constant  pressure  is 


C  -  c 
P  v 


3(3\+2G)ja  T 


1  + 


7c 


=  c  *  c 

V 


(2.15) 


As  a  result  of  these  simplifications  Eq.  (2.12)  can  be  written  simply  as 


o2o  .  b2a 

X  1  X 


3x2  a2  St2 


a 

gc  St  , 


(2.16) 


with  the  heat  equation  reduced  to  the  familiar  result 


3t 

q  -  yc  St 


(2.17) 


Hence,  Eqs.  (2.16)  and  (2.17)  are  equivalent  to  our  original  equation, 


13. 


Eq.  (2.8),  which  is  sufficient  for  problems  involving  uniform  internal 

heat  generation  in  which  the  temperature  is  automatically  known  for  any 

prescribed  heating  process.  The  above  simplifications  are  in  agreement 
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with  Boley  and  Weiner,  who  performed  a  similar  analysis  for  the  dis¬ 
placement  formulation  of  a  similar  problem.  They  considered  a  more 
general  problem  in  which  the  temperature  is  space  dependent  as  well  as 
time  dependent,  and  show  that  the  coupling  term  in  (2.9)  can  be  neglected 
if  the  dilatation  rate  is  of  the  same  order  of  magnitude  as  the  temper¬ 
ature  rate.  In  all  thermoelastic  problems  there  is  an  additional 
restriction  on  the  use  of  the  complete  energy  equation.  It  is  necessary 

that  It  -  T  I  <  T  .  However,  this  restriction  can  be  relaxed  for  the 
1  o'  o 

linear  theory  involving  small  strains  since  even  for  high  temperatures 
|T  -  Tq|  is  normally  less  than  the  strain  level  (~  10  J)  permitted  by 
the  infinitesimal  theory.  A  complete  discussion  of  this  is  given  by 
Chadwick. 

Returning  now  to  the  problem  of  the  inertial  stresses  induced  in 
a  rapidly  heated  rod,  the  equations  to  be  solved  are 


a2g  i  a2o 
ax2  '  a2  at2 


(2.8) 


(2.17) 

(2.18) 

where  the  subscript  x  has  been  dropped  for  convenience,  and 

o  •  o(x,t)  T  ■  T(t) 

u  «  u(x, t)  q  -  q(t) 


aT 

q  -  7c^ 


and 


au  o 

s  “  i  +  aT 
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In  this  work,  the  temperature  rise  will  be  kept  small  enough  so  that 
variation  in  material  properties  can  be  neglected.  However,  if  it  is 
desired  to  incorporate  temperature  dependence  of  E  and  a,  Eq.  (2.8) 
could  be  solved  by  numerical  methods.  Let  us  now  consider  the  elastic 
response  of  a  rapidly  heated  rod  for  three  special  cases  of  heating 
functions  shown  in  Fig.  2  which  have  been  chosen  for  their  similarity 
to  experimental  heating  conditions.  These  are  the  cases  of  constant 
internal  heat  generation  for  a  finite  time,  the  exponential  temperature 
rise,  and  the  exponential  temperature  rise  with  zero  initial  conditions, 
i.e.,  zero  initial  stress  rate  as  well  as  zero  initial  stress. 

(a)  Case  I  -  Constant  Internal  Heat  Generation  Rate 

Let  a  rod  be  heated  by  a  constant  internal  heat  generation  rate 
acting  for  time  t^.  This  can  be  written  as 

q  -  qQ[H(t)  -  HCt-tp]  (2.19a) 

where  H  is  the  unit  step  function,  i.e.. 


H(t-t1) 


(0  t  <  tx 

h  t  >  tx 


The  corresponding  temperature  history  is  obtained  from  (2.17),  and  is 


T  - 


m 


(H(t)  +  (i  - 1-) 


(2.19b) 


where  Tm  *  qQ/7c  ■  maximum  steady-state  temperature.  These  functions 
are  shown  in  Fig.  2a.  Since  the  stress  is  proportional  to  the  temperature, 
which  has  a  finite  rate  at  zero  time,  we  cannot  specify  a  zero  initial 
stress  rate.  Hence,  the  initial  conditions  are 
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c(x,o)  -  0 

6(x,o)  «  bQ 


(2.20) 


Also,  since  Che  heating  is  considered  to  be  rapid,  the  temperature  is 

independent  of  x.  This  implies  that  Oq  is  a  constant  and  that  the 

stress  will  reach  its  maximum  value,  -  E  a  T  ,  in  time  t, .  Hence 

m  i 


. 

a 

o 


E  a  t_ 


m 


(2.21) 


where  agt  is  the  stress  developed  if  the  bar  were  clamped  at  both  ends 

and  heated  slowly  to  T  .  This  initial  condition  will  be  verified  later 

m 

from  the  solution  of  (2.8).  In  order  to  solve  (2.8)  Laplace  transform 
methods  are  used,  and  Eq.  (2.8)  in  the  transform  o(x,s)  is 


~  -  ~  (s2?  -  sO(o)  -  o(o))  -  p  a [s2T  -  sT(o)  -  T(o)]  (2.22) 

dx  &L 

where  x  is  the  Laplace  transform  parameter,  and 


JL  Co(*»t)J  -  o(x,s)  -  /  e  o(x,t) 


dt 


T  -  T 


_L_. 

-st,  . 

1  -st, 

0  r  1 

_L.  . 

r.  '8til 

1  -  e 

2  + 
Ls  tx 

s  e 

^8tl  s2t/J 

“  Tm 

s2tj  . 

and 


T(o)  -  0 


T(o) 


Substitution  into  (2,8)  yields 

.2_  -2  0Ofc  POT  -st. 

do  s  St  _ m  1 

I?’?'  ■  h  *T  * 


(2.23) 
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The  solution  e£  this  non-homogeneous  ordinary  differential  equation  is 

-8t, 

sx/a  -sx/a  2 


o(x,s)  *  C^e  +  Cje  -  — 


°St  P°V 


’In 


a2t. 


(2.24) 


The  arbitrary  constants  and  C2  can  easily  be  evaluated  from  the 

boundary  conditions.  Since  the  bar  is  unrestrained 


a  (i,s)  -  a  (-i,s)  -  0  , 

substitution  into  (2.24)  yields 


C1  "  C2  " 


1 

r  2  '8tll 

oSt-paa  Te 

2  cosh  si/a 

-  a 

J 

(2.25) 


Since  poa  T  *  pa( E/p)T  *»  EOT  ■  a0„  ,  substitution  of  (2.25)  in  (2.24) 
m  m  m  at 

gives 


cf(x,s)  «  o, 


St 


“  - 

r  -st,‘ 

cosh  sx/a  ^ 

1  -  e  1 

cosh  si/a 

.  s2tl  J 

(2.26) 


In  order  to  invert  this  solution  back  to  the  time  domain,  we  note  the 
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following  inverse  Laplace  transforms: 

-1 


£l 


1  cosh  sx/a 


2  cosh  si/a 


-at  _ 
e  f(s) 


1  x  flat 

I 


.  U  Z  (-1)  ,  lv«x  .  ,  1  x  fl 

■  t  +  — j-  I  — — I  c°s(n  “  2)j~  sln(n  “  2  '  7 

«  a  n-1  (n  -  -  ) 

(2.27a) 

[f(t-t1)]H(t-t1)  .  (2.27b) 


Rewriting  (2.26), 


°St  h 


r  -st 

1  cosh  sx/a  _  /I  cosh  sx/a  \  cl_  1  -  e 

2  cosh  si/a  \  2  cosh  si/a  / e  2 


17. 


and  making  use  of  (2.27)  and  (2.28)  we  obtain 


2J 


St  s  at.  n-1  (n  -  -r  ) 


(-1)  .  1  .  sx  .  ,  1  * 

-  2  c°s(n-  j  )  j-  sin(n  -  ^  ) 


nut 


2i 


s2at. 


00 
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n-1  (n  -  |  ) 


(-l)n  .  1  .  nx  .  ,  lx  sa(t-t, ) 

^ — jr  cos(n  “  2  >  T”  8in^n  ~  1'  £  1 


H(t-tx) 


Defining  the  fundamental  period  of  free  vibration  of  the  bar  as  f  , 
where 
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(2.28) 


then  we  can  rewrite  this  solution  in  terms  of  the  ratio  tj/x  •  Hence, 


For  t  <  tj 


St  2s 


?{k)l 


■^frr  cos<n  -  \  >  f  Bin  “n  C’ 


(n  -  7  ) 


(2.29a) 


and  for  t  >  t^ 


£ _ L(L_)  £ 

ffSt  2n2  'C1 '  n-1 


(-1)  ,  1  x  sx 

itf*08^  *  2  >  r 


(n  -  f  ) 


sin  cn  t  -  sin  co  (t-t.k 
n  n  i  1 


(2.29b) 


where  -  (n  -  l/2)  -j  -  4s  (n  -  l/2)  ™  .  The  displacement  u(x,t)  can 
readily  be  obtained  from  an  Integration  of  Eq.  (2.18)  with  the  condition 
that  u(o,t)  -  0.  The  result  is 


St 


T"  H(t)  + 
C1 


n  -f 

C1 


t-ti>  j  +  (:H"  3  8tn(p  - 1  >  f 
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sin  ffl  t  - 
n 


(2.30) 


18. 


where  u_„  -  i  Of  T  ■  maximum  static  displacement  if  the  bar  were  heated 
slowly  and  allowed  to  expand  freely.  Equations  (2.29)  and  (2.30)  con¬ 
stitute  the  complete  solution  for  the  axial  stress  and  displacement. 

They  satisfy  the  boundary  and  initial  conditions  and  have  been  verified 
by  substitution  back  into  the  original  differential  equations  (2.2)  and 
(2.8). 

In  order  to  illustrate  the  physical  significance  of  the  solution 

(2.29)  it  is  noted  that  the  amplitude  is  solely  a  function  of  the  ratio 

tj/r  .  This  is  the  ratio  of  the  heating  time  t^  to  the  characteristic 

mechanical  response  time  t.  Hence,  it  is  possible  to  calculate 

a  /o„^  and  u  /u„^  as  a  function  of  t./t.  It  should  be  noted 
max'  S t  max'  St  1' 

that  a  cursory  glance  at  (2.29)  indicates  that  as  (tj/r)-»-0, 

(a/ogt)  -*-<*>.  This  is  not  the  case,  however,  since  a  limiting  process 
must  be  followed.  If  t^  0,  then  for  t  <  t^,  sin  (to^t)  -+-0,  and 
for  t  >  tx 

(l-cos<n  t,)sincn  t  +  sino)  t.  cosco  t 
_ n  l _ n _ n  1 _ n 


Expanding  the  sine  terms  in  their  series  form  and  passing  to  the  limit 
we  obtain 

oo  n 

«  7  2  -(-4  ,  cos(n  -  k  )  T’  coso)  t  .  (2.31) 

tj-*-0  St  ,  n«l  (n  -  j  ) 

Hence,  a  _  /o_.  -  ±  1.0  when  t./x  0.  a  _  /a-,  has  been  calculated 
for  various  values  of  tj/x  and  the  result  is  shown  in  Fig.  3.  In  order 
to  illustrate  the  stress-wave  propagation  phenomena  resulting  from  rapid 
heating,  the  solution  (2.29)  has  been  evaluated  for  the  special  case  of 


2 _ 3,  r.  -fcU 

x 


St 


2  12 
2jt  n«l  (n  -  4  J 


t  1  V  «x 
cos(n  -  j  )  j~ 
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a  0.035  x  0.500  x  10  inch  304  stainless  steel  rod  subjected  to  a  constant 
Internal  heat  generation  rate  of  65,000  BTU/sec  acting  for  a  time  duration 
of  4  microseconds. 


T 


4i  _ _ 4x5 

a  0.203  x  106 


98.5  nsec. 


T  -T 
m  o 


yo  1  (65.0001 £4  x  10  °) _ 

Ikljc  =  2(0.035) (0. 5) (5) (0.28)(0. 12) 


u„,  -  i  a  (T  -T  )  -  5(9  x  lO  6) (46)  =■  0.002  in. 

o  t  mo 

agt  -  E  a  (Tm-To)  »  (30  x  106) (9  x  10‘6)(46)  »  12,500  psi. 


Since  tj^  «  t,  amax  *  °St;  hence,  for  a  temperature  rise  of  only  46°, 
it  is  evident  that  a  fairly  large  stress  can  be  developed  provided  the 
heating  time  is  short  relative  to  the  mechanical  response  time.  The 
stress  distribution  over  the  bar  length,  -  i<x<+i,  is  shown  in 
Fig.  4  for  several  successive  times,  and  in  Fig.  5  the  stress  histories 
are  shown  at  points  x  -  0,  i/4,  i/2,  and  3i/4.  The  displacement 

history  of  the  end  of  the  rod,  x  =  /  is  shown  in  Fig.  6.  These  results 
were  calculated  from  Eqs.  (2.29)  and  (2.30). 
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(b)  Case  II  -  The  Response  of  an  Unrestrained  Rod  to  an 
Exponential  Temperature  Rise 

The  temperature  function  given  in  Case  I  is  physically  unrealistic 
since  discontinuities  at  t^  are  unlikely.  Regardless  of  how  a  body 
is  heated,  the  temperature  history  can  normally  be  approximated  by  an 
exponential  temperature  rise  of  the  form  (see  Fig.  2b) 

-t/t 

T  -  Tm  (1  -  e  *)  .  (2.32) 


20. 


Again  we  specify  the  initial  conditions  as 


<J(x,o)  ■  0 

o(x,o)  *  -  oSt/t1 


(2.33) 


Substitution  in  (2.8) 


a2q  i  a2g 

ax2  ■  a2  at2 


> 


and  with  the  above  initial  conditions  this  equation  in  the  transform 
o(x,s)  is 


Using  the  same  boundary  conditions  as  before  and  following  the  same 
procedure,  the  solution  for  o(x,s)  is 


The  only  inversion  that  poses  a  problem  is  the  second  term  on  the  right- 
hand  side.  However,  using  (2.27)  and 
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we  can  make  use  of  the  convolution  integral, 
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Noting  again  that  poa  T  «  EOT  «  a.  ,  the  inversion  is 
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where  again 
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The  displacement  can  be  calculated  in  the  same  manner  as  before,  and  is 


rtx  r 


(1  -  e 


St 


‘t/Tlx  x  2  “  (-Dnsin(a-f  )  ~t 

}J  Vn-l(»-A)2[l«V’  L 


cos  (0  t  +  a)  t,  sin  o>  t 
n  n  1  n 


■t/tx 


(2.35) 


Equations  (2.34)  and  (2.35)  represent  the  complete  solution  for  Case  II, 
and  have  been  verified  by  substitution  back  into  the  original  differential 
equations  (2.2)  and  (2.8).  The  nature  of  this  solution  is  shown  in  Figs. 

7  and  8,  where  the  stress  histories  are  plotted  at  x/4  -  0  and  x/4  -  l/2 
and  the  maximum  amplitudes  are  plotted  as  a  function  of  the  ratio  t  j/t, 
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respectively.  It  should  also  be  noted  that  as  t^  0,  Eq.  (2.34) 
agrees  with  the  corresponding  solution  (2.31)  for  Case  I  when  0. 


(c)  Case  III  -  The  Response  of  an  Unrestrained  Rod  to  an  Exponential 

Temperature  Rise  Modified  for  Zero  Initial  Stress  Rate 

Both  the  temperature  functions  chosen  for  the.  previous  two  cases  are 
perhaps  unrealistic  since  they  assume  non-zero  initial  stress  rates  or 
particle  velocities,  and  this  implies  an  infinite  initial  acceleration. 

The  exponential  temperature  history,  which  is  reasonable  for  t  >  0, 
requires  some  modification  to  incorporate  zero  initial  conditions  while 
retaining  the  exponential  character  during  the  temperature  rise.  This 
can  be  accomplished  by  the  function  (see  Fig.  2c) 

”  t/pti 

T  -  Tq  =  Tffl  [1  -  (1  +  t/ptj  )  e  3  (2.36) 

in  which  p  is  selected  arbitrarily  on  the  basis  of  the  degree  of 
approximation  to  the  simple  exponential  function  given  by  Eq.  (2.32). 

A  comparison  of  (2.36)  with  (2.32)  and  the  Case  I  temperature  rise  (2.19b) 
is  shown  in  Fig.  9  for  various  values  of  the  parameter  p.  It  is  apparent 
that  if  p  *  l/3  good  agreement  with  (2.19b)  is  obtained,  and  if 
P  •  l/2  or  3/4,  good  agreement  with  (2.32)  is  obtained  for  t  >  0. 

Assuming  Eq.  (2.36)  as  the  temperature  function,  the  initial 
conditions  are  given  by 


T(o)  *  T(o)  *  0(x,o)  *  6(x,o) 

so  that  the  transformed  Eq.  (2.8)  is 
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Following  the  same  procedure  as  before,  we  find 
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we  obtain  the  solution  for  the  stress  as 
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The  displacement  is,  therefore,  given  by 


(2.37) 


St 


Pt, 


X 

i  *  # 


2  «  (-l)nsin(n  "  )  ~ 


+  =2  Z 


c  *  .  1.2  r.  .^2  2  2 1 

n  n-1  (n  -  j)  [1+0  tj  con  ] 


-  — —  2  2  2  / 

t  Ptj  2ptjCDfl  sin  co^t  +  (1-P  t^  u>n  )  ^cos 


t_ 

Pt, 


CD 


i  ♦  *v».J 


t  - .  ‘j 

1 -  j(2. 


38) 


where,  as  shown  before, 
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The  Eqs.  (2.37)  and  (2.38)  represent  the  complete  solution  for  Case 
III  type  heating,  and  have  been  verified  by  substitution  back  into  the 
original  differential  equations  (2.2)  and  (2.8).  This  solution  also 
agrees  with  (2.31)  for  an  Instantaneous  temperature  rise,  i.e.,  t^  ■  0 
or  P  *  0.  Equation  (2,37)  is  compared  with  (2.34)  in  Fig.  10  for  the 
special  case  of  t^/r  =*  0.1,  and  p  =  3/4.  Since  the  two  are  signifi¬ 
cantly  different  only  during  the  initial  compression  pulse,  or  transient 
part  of  the  solution,  but  agree  reasonably  well  during  the  steady-state 
oscillation,  the  stress  histories  at  a  point,  as  shown  in  Fig.  7,  need 
not  be  constructed. 

The  significant  feature  of  all  of  these  solutions  is  the  amplitude 

dependence  upon  the  ratio  t ^/t.  For  Case  III  heating,  this  is  shown 

in  Fig.  11  for  P  *  l/ 2  and  P  -  3/4.  It  should  be  mentioned  that  all 

of  the  plots  of  the  theoretical  solutions  in  this  work  were  obtained  from 

calculations  performed  on  the  IBM  650  Digital  Computer  at  the  Lawrence 

Radiation  Laboratory  at  Livermore,  California,  The  series  solutions 

were  evaluated  by  calculating  as  many  terms  in  the  series  as  necessary 

-4 

to  reduce  the  partial  sums  to  magnitudes  on  the  order  of  10  ,  Hence, 

the  theoretical  curves  shown  in  this  work  represent  an  accuracy  of  about 

0.01%. 

(d)  Summary  of  Theoretical  Calculations 

For  the  purposes  of  clarity  the  following  tabular  summary  is 
prepared.  Associated  with  this  summary,  Fig.  12  shows  the  comparison 
of  the  amplitude  dependence  upon  t^/T  for  the  three  different  types 


of  heating  functions. 
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3.  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURES 

(a)  General 

Investigations  of  the  inertia  effects  of  rapid  heating  are  extremely 
difficult  to  perform  experimentally.  The  primary  difficulty  is  obtaining 
very  high  heat  generation  rates  or  thermal  fluxes  which  can  be  applied 
in  a  very  short  time  (see  Figs.  3,  7,  and  11).  Heat  sources  such  as 
solar  energy,  plasma  jets  and  propellants  produce  high  thermal  fluxes, 
but  the  heating  times  are  necessarily  long  (on  the  order  of  10  *  sec) 
demanding  excessive  natural  periods  of  the  heated  test  specimen.  Also, 
since  only  one  surface  is  heated,  it  becomes  necessary  to  determine  the 
Instantaneous  temperature  gradient  in  order  to  evaluate  the  inertial 
effects  separately.  Further,  any  theoretical  solution  of  Eq.  (2.8) 
for  which  the  temperature  function  is  also  spatial,  as  well  as  time 
dependent  poses  analytical  difficulties  which  complicate  verification 
of  experimental  results.  Additional  practical  difficulties  arise  in 
the  design  and  operation  of  high  thermal  radiation  sources. 

It  was  decided,  therefore,  to  make  use  of  electrical  heating 
methods.  Joulean  heating  will  produce  a  reasonably  uniform  temper¬ 
ature,  and  since  it  is  dependent  upon  the  current,  can  be  easily 
measured.  In  order  to  produce  short  heating  times,  however,  the 
electrical  system  must  have  a  low  Inductance  so  that  the  energy 
deposition  time  is  short.  Also,  since  very  high  heat  generation 
rates  are  required,  the  system  must  have  a  "high"  capacitance 
operating  at  a  high  voltage.  These  requirements  are  similar  to 
those  faced  in  investigations  of  exploding  wire  phenomena,  plasma 
physics,  and  magnetic  metal-forming  studies.  Researchers  in  these 
activities  at  the  Lawrence  Radiation  Laboratory  have  designed  and 
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built  several  types  of  low- inductance,  high  energy,  mobile  capacitor 
banks  which  have  an  inductance  on  the  order  of  10  ^  henries  and  an 
energy  storage  capability  as  high  as  50  kilojoules.  Depending  on  the 
external  circuit  resistance,  discharge  times  for  these  systems  are  on 
the  order  of  microseconds.  These  units  are  compact,  simple  to  operate 
and  maintain,  and  provide  reproducible  energy  deposition  rates.  Hence, 
the  loan  of  one  of  these  units  was  requested  for  use  in  this  work. 

(b)  Specimen  Geometry  and  Heating  Method 

The  experimental  problem  consisted  of  determining  the  dynamic 
elastic  strain  in  a  long,  thin  unrestrained  rod  subjected  to  rapid 
internal  heat  generation  such  that  the  heating  time  was  short  relative 
to  the  fundamental  natural  period  of  oscillation  of  the  rod.  The 
theoretical  solution  for  the  behavior  of  an  unrestrained  or  "free-free" 
rod  was  first  worked  out  as  shown  previously  in  Section  2.  However,  a 
completely  "free-free"  rod  is  dificult  to  simulate,  but  due  to  the 
spatial  symmetry  exhibited  by  the  solution  (see  Eqs.  (2.34)  and  (2.35), 
for  example),  the  results  for  a  bar  of  length  fixed  at  one  end  and  free 
at  the  other  are  identical  to  the  results  for  a  "free-free"  bar  of 
length  2.  Hence,  the  test  specimen  was  chosen  to  be  a  "fixed-free"  rod 
clamped  at  one  end  and  submerged  in  a  shallow  mercury  pool  at  the  other 
for  the  necessary  electrical  continuity.  Heating  was  accomplished  with 
a  Lawrence  Radiation  Laboratory  "Mosquitoette"  capacitor  bank  with  a 
capacity  of  14  pfd,  an  inductance  of  about  0.1  ph,  an  internal  resistance 
of  0.02  ohms,  and  with  a  voltage  control  system  for  charging  to  any 
voltage  up  to  20  kilovolts  so  that  the  maximum  energy  storage  capacity 
is  2800  joules. 
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The  heating  jig  is  shown  schematically  in  Fig.  13.  This  arrangement 

was  chosen  in  order  to  equalize  the  magnetic  forces  on  the  specimen,  and 

to  allow  control  of  the  circuit  inductance,  if  necessary,  by  varying  the 

distance  between  the  return  conductor  rods.  This  allows  some  control 

over  the  frequency  of  the  current  in  the  bar.  The  jig  can  hold  specimens 

up  to  15  inches  long  and  0.5  inches  wide. 

The  specimens  were  fabricated  from  Type  304  stainless  steel  and  for 

reasons  discussed  later  were  rectangular  in  cross-section  with  dimensions 

of  0.050  x  0.250  inches  or  0.125  x  0.250  inches  with  lengths  from  about 

2 

1.50  to  12.0  inches.  In  order  to  maximize  the  i  R  heating  and  the  time 
for  heat  conduction  out  the  ends,  stainless  steel  was  chosen  for  its 
high  electrical  resistivity  and  low  thermal  conductivity.  Also,  all 
other  physical  and  mechanical  properties  are  well  known. 

(c)  Electrical  Heating  Measurements 

The  system  described  thus  far  consisted  simply  of  an  R-L-C  circuit. 
The  current  in  the  bar  will  be  governed  by  the  circuit  constants,  and 
for  the  oscillatory  case  is  given  by 
CV  ~  .  -at 

i  *  (a  +  b^)  e  sin  bt  (3.1) 

where 

L 

_  (3.2) 

1_  R*_ 

“  4L2 

Since  the  initial  charge  voltage  VQ  and  the  capacitance  C  are 
known  (the  capacitance  of  the  jig  can  be  neglected  relative  to  14  nfd), 
it  was  necessary  to  determine  only  the  inductance  L  and  total  resistance 
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R  to  evaluate  the  current  and  hence  the  heating  time.  This  was  accom¬ 
plished  with  an  inductive  probe  consisting  of  a  small  coil  wound  around 
a  glass  tube  placed  next  to  the  specimen.  The  current  induced  in  the 
coil  was  measured  with  a  Type  551  dual-beam  Tektronix  oscilloscope  fitted 
with  Type  L  preamplifiers.  The  output  of  the  coil  is  proportional  to  the 
first  time  derivative  of  the  current  in  the  bar  (or  the  oscillating 
magnetic  field);  hence,  an  RC  integrator  was  placed  in  series  with  the 
probe.  The  non- integrated  signal  was  also  recorded  as  a  precautionary 
measure.  This  circuitry  is  shown  in  Fig.  14.  From  Eq.  (3.1)  the  period 
of  oscillation  of  the  current  is 


and  the  amplitude  decay  characteristic  can  be  obtained  from 


R_ 

2L 


*°8  ij/i2 


(3.4) 


where  ij  and  i^  are  two  successive  maximums  of  the  same  sign.  Thus, 
a  measurement  of  the  period  of  oscillation  and  the  logarithmic 
decrement  of  the  probe  signal  allows  a  calculation  of  L  and  R  from 
(3.3)  and  (3.4). 

In  order  to  calculate  the  characteristic  heating  time  it  was  assumed 
that  the  current  was  uniform  throughout  the  specimen.  Actually,  however, 
due  to  the  skin  effect  the  current  must  be  greater  at  the  surface.  This 
effect  was  calculated  from  standard  formulae,  and  the  aforementioned 
specimen  dimensions  chosen  so  that  the  effective  increase  in  resistance 
was  on  the  order  of  1.1  greater  than  the  D-C  value  (the  bar  dimensions 
are  roughly  twice  the  skin  depth).  The  temperature  history  due  to  the 
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Joulean  heating  is  then  obtained  from  Eq.  (2.17) 


7  c 


dT 

ST 


where  R  is  the  resistance  of  the  specimen  alone.  Substituting 
6 

Eq.  (3.1)  and  integrating,  the  temperature  history  is 


(3.5) 

This  result  assumes  that  no  heat  is  conducted  or  radiated  away,  which  is 
a  valid  assumption  relative  to  the  microsecond  time  scales  which  occur 
in  this  work.  This  is  discussed  in  detail  in  Sec.  2  of  the  Appendix. 

A  comparison  of  this  temperature  history  with  the  simple  exponential 
(2.32)  and  the  compound  exponential  (2.37)  is  shown  in  Fig.  15.  It  is 
evident  that  the  sine  terms  in  the  parenthesis  in  Eq.  (3.5)  are  small 
compared  to  unity  which  gives  the  approximation  (see  Case  II,  Sec.  2) 

*  -  T,  •  V1  '  *  > 


Or,  to  Incorporate  zero  initial  temperature  rate  conditions  (Case  III, 
Sec.  2)  t 


T  -  T 
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CV 


T  - 
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(-)  ^ 

\  Wc  /  R 
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9.48  x  10 BTU/ joule 


The  electrical  method  of  heating  is  obviously  very  convenient  since 
the  characteristic  heating  time  t^  can  be  controlled  by  the  inductance 
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and  resistance.  With  the  freedom  to  vary  the  length  of  the  bar,  the 
mechanical  response  time  t  can  be  varied  (see  Eq.  [2.28]),  which  then 
allows  experimentation  at  any  desired  ratio  of  tj/t.  Also  within  the 
energy  capability  of  the  capacitor  bank  and  the  heat  capacity,  Wc,  of 
the  bar,  the  maximum  temperature  can  be  closely  controlled.  However,  no 
means  of  direct  instantaneous  temperature  measurement  was  provided.  This 
is  the  major  restriction  in  this  work  and  its  significance  will  be 
discussed  later. 

(d)  Measurement  of  Axial  Strain  in  the  Heated  Rod 

The  most  direct  method  of  strain  measurement  is  to  scribe  a  grid 
on  the  specimen  surface  and  observe  the  resulting  movement  with  a  high¬ 
speed  camera.  However,  since  the  present  study  is  restricted  to  elastic 
motion,  the  strains  are  very  small  and  cannot  be  measured  directly  in 
this  manner. 

An  alternate  method  is  to  fasten  strain  gages  directly  to  the  bar, 
and  observe  their  output  electronically.  Since  electrical  discharge 
provides  the  only  convenient  method  of  heating  the  bar,  any  electronic 
measurement  or  equipment  must  be  shielded  or  removed  from  the  high 
transient  magnetic  field  surrounding  the  specimen  during  the  high  voltage 
discharge  of  the  capacitor  bank.  The  millivolt  output  of  a  stain  gage 
system  mounted  on  the  bar  could  not  be  adequately  shielded  to  obtain  an 
adequate  signal  to  noise  ratio.  However,  this  technique  was  useful  for 
determining  the  steady-state  oscillations  observed  after  the  heating 
pulse,  and  was  used  with  reasonable  success  at  low  bank  voltages.  Above 
bank  voltages  of  about  5  kv,  arcing  occurred  and  destroyed  gage  continuity. 
It  has  been  suggested  that  the  heated  bar  itself  could  be  used  as  a  strain 
gage  by  measuring  the  change  in  resistance  over  a  section.  In  addition 
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Co  Che  problem  of  shielding,  chis  method  would  measure  Che  InCegraCed 
change  in  resiscance  over  che  bar  giving  only  Che  InCegraCed  scrain. 

Since  Che  inscancaneous  currenCs  in  che  bar  can  be  as  high  as  75,000  amps, 
ic  was  apparenC  Chat  opcical  methods  were  more  likely  Co  be  successful  for 
measuring  Che  small  elastic  scrains  during,  as  well  as  afCer,  Che  heating 
pulse,  and  thac  electronic  methods  are  rendered  unsatisfactory  by  magnetic 
shielding  difficulties. 

Since  Che  elastic  strains  in  the  heated  rod  are  very  small,  some 

means  was  sought  to  amplify  their  effect.  In  che  technique  known  as 
34  35  36 

"Photostress"  *  *  the  bar  is  coated  with  a  birefringent  material 

which  produces  interference  fringes  as  the  parent  bar  is  deformed.  Using 
a  reflecting  polariscope  with  an  intense  light  source,  the  fringe  pattern 
could  be  photographed  with  a  high-speed  camera.  Since  the  coating  must 
be  necessarily  thin  to  prevent  restriction  of  the  bar  movement,  and 
since  the  fringe  order  and  the  resultant  accuracy  are  directly  propor¬ 
tional  to  the  coating  thickness,  this  method  proved  to  be  too  insensitive. 
Hence,  this  technique  was  modified  in  the  following  manner. 

A  0.25  inch  thick  by  1.0  inch  deep  beam  of  a  photoelastic  material 
(Bakelite  ERL  2774)  was  loaded  in  pure  bending  in  the  plane  of  the  beam. 
The  resulting  fringe  pattern  was  frozen  into  the  loaded  beam  by  raising 
the  temperature,  holding  it  constant  at  165°F  for  four  hours,  and  then 
allowing  it  to  cool  slowly  back  to  room  temperature.  When  viewed  through 
a  crossed  circular  polariscope  the  residual  fringe  pattern  appeared  as 
shown  in  Fig.  16.  The  beam  was  polished  on  its  lateral  surfaces  and 
then  cut  into  0.050  inch  strips  which  were  glued  to  each  side  of  the 
rectangular  specimen  as  shown  in  Fig.  17.  An  epoxy  resin  glue  commonly 
used  in  "Photostress"  techniques  was  found  satisfactory  for  this 


33 


application.  If  careful  attention  was  given  to  the  cutting  technique, 
the  fringe  pattern  remained  unaltered.  In  order  to  eliminate  thermal 
effects  during  cutting,  it  was  necessary  to  perform  all  cutting  and 
polishing  by  hand.  Polishing  of  the  surface  to  be  glued  to  the  bar  was 
done  in  successive  stages  by  using  emery  cloth  of  decreasing  grit  size 
mounted  on  plate  glass  since  it  was  important  to  maintain  a  plane 
surface  in  order  to  eliminate  initial  bending. 

With  the  photoelastic  coating  mounted  in  the  manner  shown  in  Fig. 

17d  it  was  possible  to  use  incident  rather  than  reflected  light  as  is 
normally  done  in  "Photostress"  techniques.  Also,  the  initial  fringe 
pattern  allows  accurate  determination  of  fractional  fringe  order  changes, 
and  by  observing  the  movement  of  the  entire  fringe  pattern,  a  complete 

stress  distribution  over  the  one  inch  coating  length  can  be  obtained. 

37 

This  is  discussed  in  some  detail  by  Zandman  who  has  recently  contri¬ 
buted  some  interesting  variations  of  the  photoelastic  stress  gage.  In 
order  to  observe  the  fringe  motion  during  heating,  a  Model  601  high¬ 
speed  framing  camera  was  used  with  a  crossed  circular  polariscope  and 
white  light  so  that  a  dark  field  with  integral  fringe  orders  was  obtained. 
A  high  intensity  light  source  was  designed  which  consisted  of  an  Edgerton, 
Germeschausen  and  Grier  Type  FX-1  Xenon  flash  tube  mounted  at  the  focus 
of  a  cylindrical  parabolic  reflector.  A  standard  Lawrence  Radiation 
Laboratory  (IWL-1733A)  pulse-line  chassis  was  modified  to  provide  power 
for  the  tube  as  shown  in  Fig.  18.  This  system  provided  approximately 
250  (isecs  of  sufficient  light  with  a  2  nsec  rise-time  to  expose  35  mm 
Kodak  high  speed  Ektachrome  film  with  camera  speeds  of  approximately 
350,000  frames/sec  at  f:23. 


Since  it  was  considered  sufficient  to  obtain  the  stress  at  only  one 
point  on  the  bar  (at  x  *  $/2),  white  light  was  used  which  aided  the 
identification  on  the  film  record  of  the  zero  and  two  first  order  fringes. 
The  zero  order  is  always  black  and  the  two  first  order  fringes  are  always 
marked  by  the  passage  from  red  to  blue.  However,  higher  fringe  orders  in 
the  initial  pattern  cannot  be  observed  easily  with  white  light.  If  the 
stress  distribution  over  the  entire  coating  length  is  desired,  it  is 
necessary  to  use  a  5461A®  green  filter  with  black  and  white  film  to 
clearly  mark  the  higher  order  fringes,  as  shown  for  example  by  the 
static  photograph  in  Fig.  16.  The  reasons  for  using  only  white  light 
will  be  discussed  later. 

The  experimental  apparatus  used  in  an  actual  test  run  is  shown  in 
Fig.  19.  The  light  from  an  exploded  bridge  wire  mounted  next  to  the 
specimen  marks  "zero"  time  on  the  photographic  record  and  also  by 
electromagnetic  induction  causes  a  pulse  to  be  superposed  on  the 
oscilloscope  trace  of  the  output  of  the  current  probe.  In  this  manner, 
the  photographic  record  is  synchronized  with  the  electronic  record  of 
the  current.  A  535  Tektronix  oscilloscope  fitted  with  a  Type  H  pre¬ 
amplifier  was  used  to  measure  the  camera  rotor  speed  by  recording  the 
output  of  a  reluctance  pickup  mounted  on  the  rotor  shaft.  This  accurately 
fixes  the  interframe  time  of  the  camera  record.  Hence,  a  synchronized 
record  of  the  fringe  pattern  in  the  coating  is  obtained  simultaneously 
with  the  current- time  relationship.  The  schematic  block  diagram  and 
the  function-time  diagram  for  this  procedure  are  shown  in  Figs.  20  and 
21,  respectively.  A  photograph  of  a  specimen  mounted  in  the  heating  jig 
is  shown  in  Fig.  22.  It  should  be  noted  that  to  compensate  for  film 
shrinkage  a  calibration  length  was  provided.  This  consisted  of  two 
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lines,  a  known  distance  apart,  scribed  on  a  piece  of  Polaroid  mounted 
next  to  the  coating.  This  allowed  light  to  pass  through  the  crossed 
polariscope  so  that  the  scribe  marks  were  illuminated.  Some  sample 
specimens  are  shown  in  Fig.  23. 

(e)  Operating  and  Data  Reduction  Procedure 

With  the  apparatus  arranged  as  shown  previously  in  Figs.  19  and  20, 
the  test  was  carried  out  as  follows.  Since  the  Model  601  camera  was 
used  in  the  "unsynchronized"  mode  of  operation,  no  synchronized  shuttering 
facility  was  available.  Hence,  the  room  containing  the  heating  jig  was 
darkened  and  the  camera  ran  at  the  desired  rotor  speed  with  the  mechanical 
shutter  held  open  manually.  Closure  of  the  "fire"  switch  on  the  control 
unit  dispatched  a  500  volt  fiducial  pulse  which  triggered  the  oscillo¬ 
scopes,  discharged  a  condenser  unit  to  explode  the  bridge-wire,  discharged 
the  main  capacitor  bank  through  the  specimen,  and  triggered  the  charged 
pulse-line  for  the  light  source.  Time  delay  units  were  included  in  this 
equipment,  but  were  not  needed  since  the  inherent  delays  in  the  system 
proved  satisfactory. 

The  oscilloscopes  triggered  first,  at  5  microseconds  the  light  source 
and  bridge-wire  discharged,  and  at  12  to  20  fiisec  the  bank  discharged.  Both 
the  light  source  and  the  bridge-wire  light  outputs  had  rise  times  of  about 
2  psec;  hence,  both  appeared  on  the  "zero"  frame  of  the  film  strip.  In 
order  to  prevent  rewriting  the  light  source  duration  was  kept  less  than 
the  time  for  one  revolution  of  the  camera  rotor  since  no  shuttering  was 
used.  A  250  usee  light  source  was  used  as  the  camera  had  164  frames  and 
interframe  times  from  2  to  3  ^second  were  desired.  The  oscilloscope 
traces  were  recorded  permanently  with  Dumont  Oscilloscope  Cameras  loaded 
with  Type  400  Polaroid  Film.  Sample  oscilloscope  traces  of  the  rotor 
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frequency  and  the  output  of  the  induction  probe  are  shown  in  Figs.  1A  and 
2A  of  Sec.  1,  Appendix.  The  use  of  this  data  is  also  shown  in  Sec.  1, 
Appendix,  where  the  complete  data  reduction  procedure  is  illustrated  for 
two  sample  runs. 

After  development,  the  35  mm  color  film  was  scanned  with  a  10X  Bausch 
and  Lomb  binocular  microscope  fitted  with  an  "x-y"  travelling  stage 
graduated  to  0.01  mm.  The  positions  of  the  zero  and  two  first  order 
fringes  in  the  coatings  were  measured  in  each  frame  with  respect  to  the 
fixed  calibration  scribe  marks.  Also,  the  known  distance  between  these 
marks  was  measured  to  provide  a  correction  for  film  shrinkage  from  frame 
to  frame.  Hence,  the  measurements  per  frame  consisted  of  six  fringe 
positions  and  a  calibration  length,  giving  a  total  of  seven  separate 
measurements.  The  average  movement  of  the  zero  and  two  first  order 
fringes  was  calculated  and  plotted  as  a  function  of  time  (or  frame 
number).  Sample  plots  are  shown  in  Figs.  3A  and  4A.  Since  the  move¬ 
ments  of  the  fringe  patterns  in  the  coatings  glued  on  each  side  of  the 
bar  were  averaged,  any  effect  due  to  lateral  bending  of  the  specimen  was 
eliminated  automatically,  which  insured  measurement  of  only  the  longitud¬ 
inal  motion  of  the  bar. 

It  should  be  noted  here  that  theoretically  it  is  possible  to  obtain 
the  transient  stress  distribution  over  the  length  of  the  coating.  This 
would  require  the  use  of  a  green  filter  in  the  light  path  which  would 
clearly  display  all  fringe  orders  from  zero  to,  in  this  case  as  shown 
in  Fig.  16,  six,  but  black  and  white  film  would  be  necessary,  which 
requires  some  means  of  identification  of  the  zero  order  fringe.  This  was 
attempted,  but  it  was  found  that  with  visual  microscope  scanning  it  was 
very  difficult  to  maintain  constant  reference  to  the  zero  order.  Hence, 
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color  film  was  used  since  the  zero  order  and  two  first  order  fringes  are 
always  clearly  marked  by  black  and  the  passage  from  red  to  blue,  respectively. 
However,  the  higher  order  fringes  above  the  first  are  not  clear,  but  for  the 
purpose  of  obtaining  the  stress  at  a  point,  observation  of  only  the  zero  and 
first  orders  were  sufficiently  accurate.  In  addition  to  the  identification 
problem  with  black  and  white  film,  it  became  evident  that  since  visual 
scanning  of  only  three  fringe  movements  resulted  in  seven  measurements  per 
frame,  or  a  total  of  about  500  separate  measurements  for  a  complete  film 
record,  observation  of  the  complete  fringe  pattern  would  be  extremely  time 
consuming  and  impractical.  This  difficulty  could  have  been  eliminated, 
however,  if  it  had  been  possible  to  obtain  an  automatic  film  densitometer 
which  would  plot  the  fringe  movements  on  a  strip  chart.  Even  with  this 
equipment,  the  ensuing  calculation  of  the  fringe  orders  with  time  would 
still  exist  and  would  require  voluminous  computations. 

Once  the  fringe  movements  recorded  on  the  color  film  were  plotted, 
the  change  in  fringe  order  at  a  point  was  calculated  graphically  from  the 
relative  positions  of  the  three  fringes  as  shown  in  Fig.  5A.  This  was 
done  for  as  many  points  in  time  as  desired.  With  the  change  in  fringe 
order  at  the  initial  position  of  the  zero  order  fringe  determined  as  a 
function  of  time,  the  stress-optic  law  was  used  to  calculate  the  stress 
as  a  function  of  time.  The  derivation  of  the  relationships  used  for  this 
purpose  are  given  in  the  following  section. 

(f)  Calculation  of  Stress  from  Experimental  Data 

The  equation  governing  the  relationship  between  the  observed  bire¬ 
fringence  in  the  coating  and  the  uniaxial  stress  in  the  heated  bar  is 
obtained  in  the  following  manner.  Referring  to  Fig.  1,  consider  the 
parent  material,  or  specimen,  to  be  in  a  state  of  biaxial  stress  such 


38. 


Chat  the  stress  0^  normal  to  the  surface  s  is  zero.  If  the  coating  c 
is  securely  bonded  to  the  specimen  and  is  of  substantially  lower  stiffness, 
then  the  strains  in  the  coating  and  specimen  are  equal  at  the  interface, 


i.e 
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where  the  subscripts  s  and  c  refer  to  the  specimen  and  coating, 
respectively.  For  the  biaxial  stress  condition  the  linear  thermoelastic 
stress-strain  relations  (2.3)  reduce  to 
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It  is  important  to  recognize  that  since  the  coating  remains  at  the 

initial  ambient  temperature  (see  Sec.  3,  Appendix)  and  is  securely  bonded 

during  the  period  of  measurement,  the  only  thermal  strain  present  is  that 

of  the  specimen  and  is  denoted  by  (1+v  )  Q!  T  in  both  (3.8)  and  (3.9). 

s  s 

Also,  since  the  theory  of  thermoelasticity  is  based  on  the  linear  super¬ 
position  of  thermal  and  mechanical  strains,  6  and  €  represent  only 

X  z 

the  effect  of  the  loading  which  in  this  work  is  the  inertia  of  the  specimen. 

However,  as  pointed  out  (see  Sec.  2,  p.  8)  a  uniaxial  stress  state  can 

safely  be  assumed  so  that  e  is  zero.  Therefore,  combining  Eqs.  (3.7), 

z 

(3.8)  and  (3.9)  gives 
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The  observed  birefringence  N  in  the  coating  is  related  to  the  state 
of  stress  by  the  stress-optic  law 
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where  o.  and  are  the  largest 

In  the  plane  of  observation.  Since 
Appendix,  that  “  cr^  -  0, 


and  smallest  principal  stresses  acting 
it  will  be  shown  later  in  Sec.  4, 
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Since  the  coating  is  subjected  to  the  thermal  strain  OffiT  in  the 
specimen,  Eq.  (3.12)  will  yield  the  total  stress  in  the  coating.  Hence 
in  order  to  determine  the  stress  due  to  inertia  effects  alone,  the  stress 
E  t  T(t)  must  be  subtracted. 
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Combining  (3.13)  with  (3.10)  will  give  the  inertia  stress  in  the  bar  as 
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This  can  be  non-dimensionalized  by  dividing  through  by  Es  T^,  or 
Cgt>  Dropping  the  subscripts. 
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If  the  bar  were  heated  slowly  to  T  ,  the  average  strain  would  be 
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since  the  interdependence  between  f  and  £  is 
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Poisson's  ratio  is  not  rate  dependent;  hence,  substitution  of  (3.14)  and 

(3,15)  into  (3.13)  gives 
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Since  =  0.35  and  *  0.30,  the  quantity  in  the  parentheses  is 
nearly  unity  and  can  be  dropped,  giving  the  result 
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where  T(t)  is  given  by  (2.19b),  (2.32)  or  (2.36). 

Since  the  presence  o£  a  shear  stress  increases  the  birefringence, 

Eq.  (3.17)  should  be  corrected  for  stress  gradients  in  the  x  direction. 
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In  order  to  account  for  the  effect,  of  a  shear  stress,  Zandman  has 
provided  the  correction  factor. 


_ N _ 

V  N2  +  4h2  (m/te:)2 

which  will  reduce  the  magnitude  of  c/o'g  obtained  from  (3.17).  However, 
in  this  work,  observations  of  the  fringe  patterns  indicate  that  £N/£x  is 
very  small.  This  can  also  be  verified  theoretically  for  the  positions 
chosen  for  measurement  of  the  strain.  Hence,  Eq.  (3.17)  is  used  directly. 
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In  order  to  evaluate  the  effect  of  rate  of  loading  on  the  optical 

properties  of  the  coating,  a  separate  calibration  experiment  was  performed 
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(see  Sec.  6,  Appendix).  As  pointed  out  by  Goldsmith,  the  strain-optic 
coefficient  may  vary  appreciably  with  loading  rate,  and  any  interpretation 
of  photoelastically  obtained  data  must  rely  upon  a  calibration  at  sub¬ 
stantially  the  same  rates.  This  was  done  on  specimen  No.  2  so  that  the 
calibration  was  carried  out  at  a  period  of  oscillation  of  approximately 
100  psec.  The  results  of  this  experiment  showed  that  the  optical 
properties  are  essentially  unaffected  at  this  rate.  Hence,  the  ratio 
of  the  dynamic  to  the  static  strain-optic  coefficients  in  (3.17)  can  be 
taken  as  unity  which  allows  (3.17)  to  be  written  simply  as 


T 
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Actually,  the  coating  reinforces  the  steel  bar  and  a  correction  for 
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this  must  be  determined.  Zandman,  Redner  and  Riegner  have  worked  out 
the  correction  factors  for  the  reinforcing  effects,  but  confined  their 
analysis  to  quasi-static  stress  states.  However,  for  dynamic  loading, 
the  inertia  as  well  as  the  stiffness  of  the  coating  must  be  considered. 
Hence,  the  analytical  procedure  of  Zandman,  ejt.al.  was  modified  to 
Investigate  the  reinforcing  effect  present  during  dynamic  loading.  This 
analysis  is  performed  in  Sec.  5  of  the  Appendix,  where  it  is  shown  that 
only  the  period  of  oscillation  is  affected.  The  corrected  period  T1 
is  given  by 

t*  -  t  (3.i9) 
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The  amplitude  of  the  stress  is  Increased  slightly  since  it  has  been  shown 
to  be  dependent  upon  tj/r  which  decreases  with  the  Increase  to  t'. 
Hence,  the  period  x  in  the  theoretical  solutions  (2.29),  (2.34)  or 
(2.37)  must  be  replaced  by  t'  for  comparison  of  these  solutions  to  the 
experimental  data.  It  should  be  noted  here  that  this  correction  is  not 
exact  since  its  derivation  was  based  on  several  simplifying  assumptions 
as  noted  in  the  analysis.  The  simplifications  however,  are  not  severe, 
and  for  engineering  purposes  this  correction  factor  is  sufficiently 
accurate. 

In  order  to  obtain  the  reported  data,  then,  Eq.  (3.18)  is  used.  If 
some  means  were  available  to  measure  temperature  during  the  heating,  the 
function  T(t)  could  be  evaluated  directly.  However,  no  accurate  method 
exists  and  it  is  necessary  to  assume  a  temperature  history.  It  was 
determined  that  a  simple  exponential  rise  gave  a  good  approximation  to 
Joulean  heating.  This  is  shown  by  the  curve  in  Fig.  16  which  was  con¬ 
structed  from  a  measurement  of  the  current  with  the  inductive  probe. 

The  maximum  temperature  as  calculated  from  (3.6)  was  also  verified  by 
mounting  a  thermocouple  on  a  bar  which  was  rapidly  heated  at  various 
bank  voltages.  Observation  of  the  temperature  during  the  first  60 
seconds  of  cooling  allowed  an  extrapolation  back  to  zero  by  plotting 
the  measured  temperatures  on  semi-log  graph  paper.  It  was  found  that 
(3.6)  was  in  agreement  to  within  3%  for  temperature  increases  up  to  80°F. 
Hence,  despite  the  fact  that  no  direct  means  of  temperature  measurement 
was  provided,  it  was  evident  that  an  assumed  temperature  history  in 
(3.18)  would  not  seriously  affect  the  validity  of  the  experimental 
results.  Since  the  maximum  temperature,  and  hence  Og(,  can  be  deter¬ 
mined  accurately  from  (3.6),  the  only  assumption  arises  when  selecting 
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an  appropriate  temperature  rise.  In  order  to  consider  upper  and  lower 
bounds,  Case  II  and  Case  III  types  of  heating  were  both  Investigated 
analytically.  Thus,  the  respective  theoretical  solutions  (2.34)  and 
(2.37)  were  compared  with  experimental  results  obtained  from 


Case  II: 
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These  equations  are  used  directly  to  obtain  experimental  values  of 
a/ogt  as  a  function  of  time.  The  fringe  order  N  and  NAvg  is 
obtained  from  the  film  records,  and  Is  determined  by  a  measurement 

of  L/R  from  the  probe  record.  The  results  are  compared  with  the 

r • 

appropriate  theoretical  solution  modified  by  the  correction  factor 
(3.20).  Complete  sample  calculations  are  shown  in  Sec.  1  of  the 
Appendix  for  two  representative  tests. 
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4.  EXPERIMENTAL  RESULTS 

The  experimental  results  were  obtained  from  tests  on  four  bars  of 
differing  sizes  with  two  of  the  bars  each  having  two  different  coating 
sizes.  All  bars  were  fabricated  from  the  same  stock  of  304  stainless 
steel,  and  all  coatings  were  cut  from  the  same  piece  of  Bakelite 
(ERL-2774).  Coating  sizes  were  varied  to  evaluate  the  effect  of  the 
coating  on  the  response  of  the  bar.  This  was  accomplished  by  selecting 
one  specimen  with  large  bar  volume  and  the  other  with  small  bar  volume 
relative  to  coating  volume.  All  specimens  were  tested  under  "fixed-free 
end  conditions  with  the  exception  of  one  specimen  tested  under  an 
approximate  "free-free"  condition  for  the  purpose  of  illustrating  the 
effect  of  imperfect  clamping.  The  bars  were  rapidly  heated  to  cause 
temperature  changes  of  approximately  30°F  for  different  bar  lengths 
selected  to  provide  a  range  of  heating  ratios,  tj/r,  varying  from 
0.05  to  0.28.  All  test  parameters  are  listed  in  Table  II,  and  the 
measurements  of  the  inertial  stresses  as  a  function  of  time  are  plotted 
and  compared  with  the  appropriate  theoretical  solutions  in  Figs.  24  to 
31.  A  comparison  of  these  results  to  the  theoretical  prediction  for 
the  dependence  of  the  initial  compressive  stress  on  the  heating  ratio 
tj/r  is  shown  in  Fig.  31.  Sample  calculations  and  data  are  given  in 
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Sec.  1  of  the  Appendix. 


TABLE  II  -  TEST  PARAMETERS 
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5.  DISCUSSION  OF  RESULTS 


Theoretical  Results 

The  interpretation  of  the  theoretical  solutions  obtained  in  Sec.  2 

leads  to  some  interesting  conclusions  concerning  the  physical  nature  of 

thermoelastic  dynamics.  The  steady  state  portion  of  the  solutions  (when 

t  »  t^)  coincide  with  the  solution  for  the  response  of  a  thin  rod 
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subjected  to  initial  compression  and  sudden  release.  Hence,  thermo- 
elastic  dynamic  problems  are  exactly  analagous  to  mechanical  vibration 
or  shock  problems.  If  an  elastic  body  is  heated  rapidly  enough,  the 
inertia  prevents  rapid  thermal  expansion  such  that  a  large  inertial 
compressive  stress  state  exists  equal  to  that  which  would  result  if 
the  body  were  restrained  and  heated  slowly.  Subsequent  unloading  takes 
place  as  a  free  periodic  oscillation  about  the  initial  stress  state  as 
shown  in  Figs.  4,  5  and  7,  for  example.  It  is  evident,  then,  that  the 
maximum  stress  possible  will  be  equal  to  the  static  or  initial  value, 
and  the  corresponding  maximum  displacement  is  twice  the  static  value. 
These  conditions  occur  only  if  the  body  is  heated  instantaneously, 
i.e.,  t^/T  “  0*  which  represents  the  extreme  conditions  of  heating. 

As  shown  in  Figs.  3,  8  and  11,  the  stress  amplitude  depends  only 
on  the  ratio  of  a  characteristic  heating  time  to  a  characteristic 
mechanical  time  which  is  the  fundamental  period  of  free  oscillation 
of  the  rod.  Hence,  the  meaning  of  "rapid  heating"  is  made  clear  since 
the  inertial  effect  becomes  important  only  when  the  heating  time  is 
significantly  less  than  the  mechanical  response  time.  It  is  also 
evident  from  Fig.  12  that  variation  in  the  imposed  heating  function 
influences  the  stress  amplitude.  This  is  expected,  since  the  time  to 
reach  a  given  temperature  will  be  shortest  for  a  constant  heating 
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rate  (Case  I)  as  compared  to  an  exponential  heating  rate  (Case  II  or  III). 

It  should  be  noted  that  these  curves  are  plotted  for  x/l  -  0.5.  Due  to 
the  spatial  distribution  of  stress,  the  maximum  amplitude  is  largest  at 
x/l  ■  0;  hence  these  curves  do  not  represent  the  maximum  stress  in  the 
entire  rod,  but  only  the  maximum  stress  at  the  position  indicated. 

In  a  real  material  there  will  be  some  damping  present  which  will 
reduce  the  subsequent  oscillations  to  zero.  The  theoretical  analysis 
presented  here  does  not  include  this  effect;  however,  the  solutions  given 
do  represent  the  maximum  stress  conditions.  If  damping  were  considered, 
its  effect  would  probably  be  manifested  by  some  type  of  exponential  decay 
function  multiplying  the  series  solutions  given  in  Sec.  2.  Since  the 
mechanical  and  physical  properties  of  304  stainless  steel  are  nearly 
constant  for  small  temperature  changes,  this  work  was  performed  by 
restricting  the  maximum  temperature  rise  to  less  than  40°F.  For  large 
temperature  changes  the  temperature  dependence  of  the  elastic  modulus, 
specific  heat  and  coefficient  of  thermal  expansion  can  be  considered  in 
the  solution  of  Eq.  (2.8),  but  would  require  numerical  solution  with 
machine  computation.  However,  if  the  temperature  rise  is  large  enough, 
then  ogt  will  exceed  the  yield  stress  and  plastic  flow  may  occur  during 
the  first  compressive  pulse.  If  this  occurs  then  the  solution  is  compli¬ 
cated  by  the  deviation  from  elastic  behavior.  However,  it  was  shown  « 

(Fig.  4)  that  if  no  temperature  gradients  exist,  all  points  are  raised 
to  some  stress  level  at  the  same  time  during  the  temperature  rise.  Hence, 
yielding  would  occur  only  during  the  initial  compression,  while  subsequent 
free  oscillation  would  occur  elastically  after  some  permanent  deformation 
had  taken  place  upon  the  first  stress  reversal. 
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The  physical  nature  of  thermoelastic  dynamics  has  interesting  experi¬ 
mental  consequences  as  it  would  be  possible  to  simulate  the  effects  of 
rapid  heating  by  mechanical  rather  than  thermal  loading.  For  instance, 
a  long  rod  could  be  heated  slowly  to  the  desired  temperature  and  then 
compressed  until  the  stress  level  reached  the  value  associated  with  the 
desired  value  of  tj/T.  Sudden  release  of  the  rod  would  result  in  free 
oscillation.  An  example  of  a  utilization  of  this  idea  is  the  calibration 
procedure  used  in  this  work  and  described  fully  in  Sec.  6  of  the  Appendix. 

This  technique  could  be  useful  for  testing  at  high  temperatures.  However, 
the  photoelastic  coating  could  not  be  used,  but  instead  high-temperature 
strain  gages  would  be  utilized  for  dynamic  strain  measurement.  Obviously, 
a  rod  is  of  little  practical  interest,  but  this  technique  would  prove 
useful  for  studies  of  the  motion  of  thin  shells  at  high  temperatures  since 
the  geometry  of  such  structures  restricts  attempts  to  experimentally 
induce  rapid  uniform  internal  heat  generation  by  conventional  methods. 

The  elastic-plastic  response  of  thin  spherical  shells  has  been  investi¬ 
gated  by  the  author,  and  it  was  found  that,  as  expected,  the  results  are 
analogous  to  the  thin  rod.  Hence,  the  mechanical  method  is  feasible  and 
would  require  heating  a  shell  slowly,  compressing  it  to  an  initial  stress 
and  observing  the  subsequent  oscillation  upon  release. 

4 

Experimental  Results 

Generally,  agreement  between  the  experimental  data  and  theoretical 
predictions  is  comparatively  good  over  the  first  half  cycle  of  response, 
but  gradually  decreases  with  time.  Since  this  subsequent  disagreement 
appears  primarily  as  a  deviation  in  amplitude  but  not  frequency,  it  was 
apparent  that  some  damping  effect  was  present.  Calibration  of  the  coating 
(Sec.  6,  Appendix)  established  that  the  optical  properties  were  essentially 


rate  independent,  and  because  the  coating  temperature  remained  constant 
during  the  period  of  observation  (Sec.  3,  Appendix),  accurate  measure¬ 
ments  of  the  dynamic  strain  were  obtained.  Since  the  internal  damping 
of  304  stainless  steel  would  not  be  sufficient  to  significantly  affect 
the  first  two  or  three  cycles,  it  was  reasoned  that  either  the  forced 
oscillation  of  the  coating  introduced  sufficient  energy  loss  per  cycle 
or  that  the  damping  was  caused  by  imperfect  clamping  at  the  fixed  end. 
Experimental  investigation  of  these  effects  revealed  that  significant 
damping  was  present  and  was  primarily  caused  by  the  clamp.  Hence,  it  is 
concluded  that  the  bar  is  not  behaving  exactly  as  predicted  by  the 
theoretical  solution  which  assumes  zero  damping,  and  that  the  photoelastic 
data  provides  accurate  measurement  of  the  response  of  the  bar.  The  experi¬ 
mental  verification  of  these  conclusions  is  discussed  in  the  following. 

The  theoretical  solutions  were  based  on  idealized  boundary  conditions, 
l.e.,  the  displacement  at  the  fixed  end  of  the  rod  was  assumed  to  be 
Identically  zero,  and  that  the  free  end  is  stress  free.  It  is  unrealistic 
to  expect  a  clamp  to  be  efficient  enough  to  satisfy  the  assumed  condition 
of  zero  rod  displacement.  If  there  is  relative  movement  between  the  rod 
and  the  clamp,  then  certainly  structural  damping  of  the  oscillations  will 
occur.  The  clamp  used  in  actual  tests  consisted  of  a  0.50  x  0,63  x  0.75 
inch  copper  block  used  to  clamp  the  specimen  to  the  upper  electrode  with 
two  4-40  steel  machine  bolts,  as  shown  in  Fig.  13.  It  will  be  demonstrated 
that  depending  upon  the  tension  in  these  bolts,  the  response  of  the  bar  can 
vary  widely.  This  was  investigated  experimentally  by  clamping  a  specimen 
in  the  heating  jig  and  preloading  it  with  a  weight  suspended  by  a  thin  wire 
fastened  at  the  free  end  of  the  specimen.  Sudden  release  of  the  load  was 
achieved  by  exploding  a  section  of  the  wire  with  the  capacitor  bank.  The 
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resulting  oscillations  were  observed  with  strain  gages  mounted  on  the  bar, 
as  shown  in  Fig.  23.  This  procedure  is  identical  to  the  calibration  pro¬ 
cedure  discussed  in  Sec.  6,  Appendix  and  illustrated  in  Fig.  7k.  The  strain 
records  so  obtained  are  shown  in  Figs.  32a,  32b  and  32c.  Test  (a)  was  run 
with  the  clamping  bolts  as  tight  as  possible,  test  (b)  with  the  bolts  under 
"medium"  tension,  and  test  (c)  with  the  bolts  under  minimum  tension,  or 
Just  "snug".  These  records  were  reproducible  and  many  more  were  obtained 
for  various  combinations  of  bolt  tension  in  the  top,  bottom,  or  both  bolts. 

The  results  show  that  as  long  as  the  bolts  are  reasonably  tight  the 
strain  records  were  all  nearly  identical,  but  for  minimum  bolt  tension  the 
oscillations  were  completely  damped  as  shown  in  Fig.  32c.  Even  for  maximum 
bolt  tension  (Fig.  32a)  damping  was  present,  as  indicated  by  the  first  two 
cycles.  Apparently,  the  oscillations  after  these  two  cycles  are  drastically 
affected  by  subsequent  erratic  motion  of  the  wire  upon  unloading,  and  are, 
therefore,  not  representative  of  the  "fixed-free"  condition.  Using  the 
first  two  cycles  of  test  (b),  a  damping  factor,  or  logarithmic  decrement, 
was  calculated  and  applied  to  the  theoretical  solution  used  for  Run  2. 

The  envelope  of  this  exponential  decay  is  shown  by  the  dotted  lines  super¬ 
posed  over  the  plots  in  Fig.  25a.  The  excellent  agreement  Indicates  that 
the  apparent  discrepancy  is  not  a  result  of  the  coating  method,  but  instead 
is  due  to  imperfect  clamping.  Hence,  the  theoretical  solutions  derived  in 
Sec.  2  should  be  modified  by  the  inclusion  of  a  damping  factor  for  purposes 
of  comparison  with  experimental  data.  Since  only  one  specimen  was  large 
enough  to  be  fitted  with  strain  gages,  the  experimental  damping  correction 
was  applied  only  to  the  corresponding  heating  test,  namely  Run  2.  The 
results  of  this  correction  are  shown  in  Fig.  25b,  which  illustrates  the 
Improved  agreement  with  the  experimental  data. 
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Further  examination  of  the  strain  records  in  Fig.  32  shows  that  the 
first  two  cycles  do  have  the  same  general  wave  form  as  predicted  theoreti¬ 
cally.  However,  it  was  noted  that  in  some  tests  the  peaks  of  the 
oscillations  were  cut  off  sharply,  apparently  by  some  action  of  the  clamp. 
This  effect  is  particularly  pronounced  for  Run  1  in  Fig.  24,  and  serves 
to  explain  the  large  disagreement  between  experimental  data  and  non-damped 
theoretical  predictions.  Also,  since  the  specimen  used  for  Run  1  was 
physically  the  largest,  its  erratic  behavior  is  expected  due  to  the  decrease 
of  clamp  efficiency  for  the  larger  specimens. 

Some  thought  was  given  to  the  redesign  of  the  clamp  in  order  to  reduce 
the  damping.  An  exploratory  experiment  was  performed  by  clamping  the 
specimen  between  two  large  steel  plates  with  a  heavy  C-clamp.  It  was 
expected  that  excessive  clamp  mass  coupled  with  extreme  tightness  would 
approximate  the  best  clamp  that  could  be  designed.  The  resulting  strain 
records  obtained  as  above  indicated  that  the  damping  was  reduced,  but 
nevertheless,  still  present.  The  logarithmic  decrement  was  reduced  from 
0.0045  to  about  0.003.  Also,  the  wave  form  was  slightly  improved  for  the 
first  cycle  with  the  peak  appearing  as  predicted,  but  the  peak  of  the 
second  cycle  was  removed.  Hence,  it  was  decided  that  complete  redesign 
of  the  clamp  was  not  Justified. 

In  order  to  further  verify  the  effect  of  clamping,  a  separate  heating 
test  was  performed  on  a  "free-free"  bar.  An  attempt  was  made,  with  some 
difficulty,  to  design  end  supports  that  would  not  restrict  specimen  motion. 
This  was  done  by  supporting  the  ends  with  styrofoam  cushions  mounted  in 
5/16  inch  grooves  cut  into  suitable  copper  electrodes.  Mercury  was  used 
to  achieve  electrical  continuity  between  the  specimen  and  the  electrodes, 
and  the  entire  assembly  was  mounted  horizontally  in  the  heating  jig  which 
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had  been  rotated  90°  in  the  plane  of  the  return  conductors.  The  specimen 
was  cut  from  the  12  inch  bar  (Run  1)  such  that  the  coating  was  at  the 
center  (x/l  =  0)  of  a  10  inch  "free-free"  rod  which  theoretically  had  an 
oscillatory  period  of  100  microseconds.  The  results  of  this  test  are 
shown  in  Fig.  29,  and  much  better  agreement  with  theory  was  obtained. 

The  experimental  data  indicates  that,  for  at  least  the  duration  of 
observation,  no  damping  was  present.  As  contrasted  with  Fig.  25,  the 
difference  in  wave  shape  is  due  to  the  difference  in  location  of  strain 
measurement;  i.e.,  x/l  -  0  instead  of  x/l  =  0.5.  This  is  also  shown 
by  Fig.  7. 

The  results  from  this  test  are  compared  with  theory  for  the  two  types 
of  heating  in  Fig.  29  (Case  II)  and  Fig.  30  (Case  III,  0  -  3/4).  A  com¬ 
parison  of  these  two  figures  shows  that  best  agreement  is  obtained  for 
Case  II  heating.  This  indicates  that  the  simple  exponential  provides 
the  best  approximation  to  the  temperature  history,  although  the  actual 
time  constant  t^  may  be  longer  than  measured  by  the  inductive  probe 
method.  This  is  likely  since  the  determination  of  L/R  was  made  by 
taking  the  average  logarithmic  decrement  of  the  current  trace  which  did 
not  always  decay  in  an  exact  exponential  manner  but  was,  nevertheless, 
sufficiently  close  to  justify  this  procedure.  For  rapid  heating  by 
capacitor  discharge,  then,  it  is  evident  that  the  simple  exponential 
is  the  best  assumption  for  the  temperature  history  since  it  will  be 
slightly  conservative  and  will  give  reasonable  accuracy  with  a  minimum 
of  calculations. 

The  only  other  source  of  specimen  damping  was  the  damping  induced  by 
the  coating  motion.  In  order  to  evaluate  this  effect,  specimens  of  the 
same  size,  but  with  different  coating  size  were  tested.  Also,  a  separate 
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calibration  test  was  performed  in  which  a  bar  was  loaded  by  a  weight,  as 
mentioned  previously,  once  with  a  coating  and  once  without  a  coating,  but 
with  strain  gages  mounted  on  the  bar.  The  results  of  these  tests  showed 
that  no  significant  difference  In  response  occurred.  Hence,  it  was  con- 
c luded  that  the  coatings  used  in  this  work  introduced  no  observable 
damping. 

In  view  of  these  additional  experiments,  and  by  virtue  of  the  excellent 
agreement  obtained  in  the  calibration  test  shown  in  Fig.  8A,  it  is  evident 
that  the  photoelastic  coating  technique  can  provide  accurate  measurement  of 
the  dynamic  strain  distribution  in  rapidly  heated  specimens.  It  should  be 
noted,  however,  that  this  technique  is  applicable  only  for  short  times  since 
all  measurements  must  be  made  before  heat  transfer  to  the  coating  can  take 
place.  This  is  not  a  serious  restriction,  however,  as  it  was  shown  in  Sec. 

3,  Appendix  that  the  temperature  of  the  coating  at  a  distance  of  only  0.005 
inches  from  the  heated  bar  remains  unaffected  for  the  first  300  micro¬ 
seconds,  the  coating  is  essentially  at  room  temperature  during  this  period. 

In  this  work  only  the  stress  at  a  point  was  measured,  but  this  technique 
is  directly  applicable  to  measurement  of  the  stress  distribution  over  the 
coating.  This  can  be  done  by  using  monochromatic  light  and  black  and  white 
film,  and/or  by  spacing  coatings  all  along  the  bar.  The  only  major 
difficulty  arises  in  carrying  out  the  multitude  of  measurements  and 
calculations  involved  in  using  more  than  three  fringes  in  each  coating. 

An  example  of  the  type  of  data  obtained  from  use  of  5461  A°  light  and 
black  and  white  film  is  shown  in  Fig.  3A  of  Sec.  1,  Appendix.  This, 
however,  is  not  a  serious  drawback  provided  such  a  time-consuming  task 
can  be  justified. 
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Since  one  of  the  primary  concerns  was  to  determine  the  validity  of 
the  theoretical  relationship  for  the  dependence  of  the  inertial  stress 
on  the  heating  ratio  t^/r,  several  tests  were  conducted  on  bars  of 
shorter  lengths.  These  lengths  were  selected  to  provide  some  variation 
in  the  heating  ratio  -  from  0.05  to  about  0.30.  From  the  results  of 
these  additional  tests  shown  in  Figs.  26,  27  and  28,  the  maximum  initial 
compression  stress  was  obtained.  This  data  is  plotted  in  Fig.  31  and 
compared  with  the  theoretical  prediction  for  Case  II  heating.  Due  to 
the  damping  effect  of  the  clamp,  these  results  should  logically  lie  below 
the  predicted  values.  The  only  large  disagreement  occurs  at  t^/t  ■  0.05, 
and  is  probably  due  to  the  erratic  behavior  peculiar  to  the  larger  speci¬ 
mens  as  shown  by  Fig.  24.  These  results,  then,  illustrate  the  validity 
of  the  theoretical  solution,  and  from  these  tests  it  can  be  concluded 
that  similar  calculations  for  bodies  of  other  geometries  will  provide 
accurate  results  for  engineering  design  purposes. 

The  primary  restriction  in  this  work  is  that  the  measurement  of  the 
transient  response  during  the  heating  actually  depends  upon  an  independent 
measurement  of  the  temperature  rise.  However,  to  this  author's  knowledge, 
no  successful  temperature  measurement  device  with  microsecond  response 
time  has  been  developed  for  use  in  the  presence  of  large  magnetic  fields. 
Since  it  was  shown  that  the  temperature  rise  resulting  from  capacitor 
discharge  is  closely  approximated  by  a  simple  exponential  function,  lack 
of  a  separate  temperature  measurement  is  not  a  serious  restriction.  Hence, 
the  determination  of  stress  from  Eq.  (3.21)  will  provide  accurate  results. 
In  fact,  a  significant  result  of  this  work  is  that  the  photoelastic  coating 
technique  can  be  used  directly  to  determine  the  temperature  by  measuring 
the  fringe  order  with  time  and  using  the  strain-optic  law.  However,  the 


55. 


temperature  dependence  of  the  coefficient  of  thermal  expansion  of  the 
heated  body  and  the  dynamic  strain-optic  coefficient  of  the  coating  must 
be  known.  This  will  require  an  extended  calibration  of  the  coating  at 
the  frequencies  expected  during  the  measurement.  The  calibration  procedure 
presented  in  this  work  provides  an  excellent  technique  since  it  can  be 
performed  at  any  desired  frequency  by  simply  varying  the  length  of  the 
specimen. 

Additional  experimentation  was  performed  in  which  the  coated  specimens 
were  lapidly  heated  to  as  high  a  temperature  as  possible  with  the  2800  joule 
capacitor  bank.  These  tests  indicated  that  for  temperatures  beyond  350°F, 
the  coatings  did  not  remain  bonded.  No  actual  photographs  were  taken  to 
determine  the  fringe  movement  at  these  high  temperatures,  but  examination 
of  the  fringe  pattern  and  glue  joint  showed  no  appreciable  change.  Hence, 
it  is  concluded  that  the  coating  technique  should  be  useful  for  transient 
measurements  at  elevated  temperatures  at  high  as  350°F. 


H 


56. 


6.  CONCLUSIONS 

The  results  of  this  study  lead  to  the  following  conclusions: 

1.  The  amplitude  of  thermally  induced  stresses  in  an  unrestrained  elastic 
rod  are  dependent  upon  the  ratio  of  the  heating  time  to  the  mechanical 
period  of  oscillation.  For  a  given  total  heat  input,  the  maximum 
stress  amplitude  results  from  a  constant  heating  rate  as  compared  to 
an  exponential  type  of  input. 

2.  It  has  been  demonstrated  that  rapid  uniform  internal  heat  generation 
can  be  accomplished  by  the  discharge  of  a  low- inductance  capacitor 
bank,  and  that  this  method  provides  high  heating  rates  with  close 
control  over  the  ratio  of  heating  time  to  mechanical  response  time. 

3.  An  experimental  technique  utilizing  a  photoelastic  coating  has  been 
developed  for  observing  the  dynamic  stress  distribution  in  a  rapidly 
heated  rod.  Since  this  is  basically  an  optical  technique,  measure¬ 
ments  made  in  short  time  intervals  are  independent  of  environmental 
conditions,  particularly  magnetic  fields,  which  normally  restrict 
the  use  of  electronic  methods.  Hence,  this  method  can  be  used  to 
obtain  accurate  measurements  of  the  transient  as  well  as  the  steady 
state  mechanical  response. 

4.  Experimental  measurements  verify  the  results  of  the  analysis  which 
was  based  on  the  uncoupled  thermoelastic  field  equations,  and  also 
indicate  that  the  heating  rate  produced  by  capacitor  discharge 
through  an  underdamped  circuit  can  be  satisfactorily  approximated 

by  a  simple  exponential  function.  In  fact,  the  photoelastic  coating 
technique  can  be  used  to  measure  transient  temperatures  provided  the 
temperature  dependence  of  the  thermal  coefficient  of  expansion  of 
the  parent  material  is  known. 
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VARIATION  OF  STRESS  AMPLITUDES  WITH  RATIO 
OF  HEATK3  PERIOD  TO  FUNDAMENTAL  PERIOD 
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FIG.  4.  STRESS  DISTRIBUTION  SEQUENCE  IN  A  RAPIDLY  HEATED  ROD  (CASE  I) 


FIG.  5.  STRESS  HISTORIES  IN  A  RAPIDLY  HEATED  ROD  (CASE  II) 
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FIG.  8.  VARIATION  OF  STRESS  AMPLITUDES  WITH  RATIO 
OF  HEATING  PERIOD  TO  FUNDAMENTAL  PERIOD 
(CASE  II) 
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FIG.  12.  THE  EFFECT  OF  TYPE  OF  HEATING  ON  STRESS  AMPLITUDE  VARIATION 
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FIG.  13.  HEATING  JIG 
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FIG.  15.  COMPARISON  OF  EXACT  TEMPERATURE  VARIATION  WITH  CASE  II 
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FIG*  16*  FRINGE  PATTERN  IN  PRESTRESSED  BEAM 


4 


HID  7971 


i 


HYD  7972 


FIG.  17.  SPECIMEN  PREPARATION 


FIG.  18.  LIGHT  SOURCE  CIRCUITRY 
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21.  FUNCTION  DIAGRAM  OF  SEQUENCE  OF 


FIG.  22.  HEATING  JIG  ASSEMBLY 


FIG.  23.  SPECIMENS 


Runl.  304  s 
0.125 


Run  2  304  st.  steel  bar  ( fixed-free) 

0.125x0.250x5.00" 


Run  3  304  st.  steel  bar  ( fixed  -  free) 
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FIG.  26.  STRESS  HISTORY  AT  A  POINT 


STRESS  HISTORY  AT  A  POINT 


FIG.  28.  STRESS  HISTORY  AT  A  POINT 


Run  6 


a. 

O 

O 

IO 

<0  * 


«  0 

CM  • 


HYD  7986 


FIG.  30.  STRESS  HISTORY  AT  A  POINT 
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la)  MAXIMUM  BOLT  TENSION] 


(b)  MEDIUM  BOLT  TENSION  | 


(c)  MINIMUM  BOLT  TENSION 


FIG.  32.  STRAIN  GAGE  SIGNALS 

HYD  7988 


APPENDIX 
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1.  Sample  Data  and  Calculations 

Teats  on  two  different  specimens  were  selected  as  representative. 
Hence  the  data  reduction  on  a  "fixed  free"  test  (Run  2)  and  the  "free* 
free"  test  (Run  6)  are  chosen  for  inclusion  here.  Figs.  1A  and  2A 
illustrate  the  oscilloscope  records  of  the  output  of  the  inductive 
probe  and  the  camera  rotor  signal,  respectively.  The  fringe  move¬ 
ments  are  plotted,  as  shown  in  Figs.  3A  and  4A.  The  fringe  order 
at  any  given  time  is  then  calculated  graphically  by  the  construction 
shown  in  Fig.  5A.  Once  the  fringe  order  is  known,  Eq.  (3.21)  is 
used  to  calculate  The  calculations  of  parameters  are  shown 

in  the  following. 

Run  2 

The  camera  interframe  time  is  obtained  from  the  oscilloscope  record  in 
Fig.  2A  where  it  is  shown  that  the  period  of  oscillation  is  500  psecs. 
Since  the  rotor  speed  is  half  the  signal  frequency,  the  rotor  turns  at 
1000  rps.  From  camera  geometry,  the  interframe  time  is  2.95  Msecs. 
From  Eq.  (3.4)  and  the  inductive  probe  record 

R/2L  -  (l/Tc)log(i1/l2)  -  (106/I3.1)l°g<2.015)  -  0.0534  x  106 

From  Eq.  (3.3) 

l  -  0.275  uh.  •,  R  -  2(0.275) (0.0534)  -  0.0294 
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For  304  stainless  steel, 

r  .  t)i/A  .  (31  x  10‘6)(5)/(0.125)(0.250)  -  0.005 

8 


Hence,  from  (3.6), 


(CV2/2)(J/Wc)(R  /R)  -  T 

o  ui 

«  (14  x  10~6) (10. 3)  x  10^) (9.48  x  10“4)(0.005)  „ 

2(0. 28) (5) (0. 125) (0.250) (0. 12) (0. 0294)  * 

0-  ■  EOT  -  (30  x  106) (9  x  106) (22.6)  -  6100  psi. 

b  t  m 

tx  *  L/R  -  10‘6/(2) (0.0534)  *  9.36  |isecs. 

T  ■  4i/a  »■  4(5)/0.20  -  100  pisecs. 


m 


1  +  (Ae/As)(Pc/P.) 
I  - 


1  +  (0.06/0. 125) (0.045/0. 28) 
1  -  (0.06/0. 125)(0.5/30) 


1.088 


x'  -  yfa  r  m  •Jl.OSB  (100) 

and 

^/t'  -  0.09 


From  the  fringe  data  in  Fig.  4A,  N 


a/a 


St 


Avg 

N  -t/9.36 

ds  “  <l  '  *  > 


105  (isecs. 


0.55. 


Hence,  a  calculation  of  N  as  a  function  of  time  will  give  a/°§t  as 


a  function  of  time. 
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Run  6 

Interframe  time  2.95  Msecs. 

R/2L  -  (106/l6.7)log(2.12)  -  0.0451  x  106 

L  •  0.50  uh.  ;  R  -  2(0. 5)(0. 0451)  -  I 

R  -  (31  x  10_6)(10)/(0. 125)(0.25)  -  0.010 

s 

T  (14)(12.56)2(9.48  x  10"S(0.010) 

m  2(0.28) (10) (0. 125) (0. 250) (0. 12) (0.045) 


m 


»/«st 


EOT  -  (30) (9) (22.1)  - 

in 


L/R  -  11.1  usees.  t 

1.10,  Tr  ■  -Jm  r  ■ 


N 

0.55 


-t/11.1 
(1  -  e  ) 


6300  psi. 

■  ZJt/k  m 


.0451 


-  22. 1°F 


100  usees. 

ti/T'  - 


108  usees.  ; 


0,016 
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RUN  2 
Current 
0.5  volt/cm. 

Voltage 
2  vo  Its /cm. 


RUN  6 

Current 
0.5  volt/cm.  — 

Voltage 
2  volts /cm. — 


-Bridge  wire  Fiducial 


FIG.  1A.  OUTPUT  OF  THE  INDUCTIVE  PROBE 
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RUN  2 


RUN  6 


FIG.  2A.  CAMERA  ROTOR  SIGNAL 
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FIG.  4A.  OBSERVED  FRINGE  T 


0  Q2  04  0.6  08  1.0 
Fringe  order,  N 


FIG.  3A.  GRAPHICAL  PROCEDURE  FOR  CALIBRATION  OF  FRINGE  ORDER 
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2.  Temperature  Distribution  In  the  Heated  Specimen 

In  order  to  determine  the  time  for  significant  temperature  drop  due 
to  conduction  out  the  ends  of  the  bar,  consider  a  rod  of  length  i  and 
thermal  diffusivity  D  with  lateral  surfaces  insulated,  and  the  ends 
held  at  zero  temperature.  If  the  bar  is  initially  at  a  uniform  temper¬ 
ature  T  ,  the  temperature  in  the  bar  for  t  >  0  is  given  by* 

h2*^ 


T(x, t) 


2 

n**  1|  3j  3|  •  •  • 


a) 


if  x  ■  i/2,  then  the  temperature  at  the  midpoint  is 


til 


T  4  _ 

—  «■  —  2 

T  Jt  toe  B 

m  f  n*  « •  t 


n*l 

2  i 

—  e 


twt^t 
2 


(2) 


The  properties  of  304  stainless  steel  are 


k  -  9.4  BTU/hr/ft2/#F/ft 

w  -  485  lb/ft3 
c  -  0.12  BTU/lb°F 

D  *  k/wc  *  0.16  ft2/hr  =■  44.5  ft2/nsec. 

For  the  shortest  specimen,  let  i  ■  1.5  inches,  Hence 


m 


4 

-  2 

n  n-1,3,5... 


til 


n-1 

2 


-4n2t 


n 


where  t  is  in  units  of  seconds. 


(3) 


*R,  V.  Churchill,  Operational  Mathematics.  McGraw-Hill,  1958,  2nd  Ed., 

p.  218. 
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The  results  of  this  calculation  give  the  temperature  history  at  the 
center  of  a  1.5  inch  bar  as 

t(secs.)  0.010  0.100  0.500  1.000 

~  0.98  0.840  0.20  0.02 

Since  all  measurements  are  made  within  the  first  300  |isecs.  after 

bank  discharge,  it  can  be  concluded  that  the  temperature  of  the  bar 

remains  constant  at  T  .  This  result  was  calculated  for  the  shortest 

m 

possible  specimen;  hence,  the  time  for  significant  temperature  drop 
will  generally  be  greater  than  0.50  secs. 
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3.  The  Temperature  History  of  the  Photoelastic  Coating 

Consider  a  thin  slab,  infinite  in  extent,  of  the  photoelastic 
material  of  thickness  J  subjected  to  a  temperature  T  maintained 

ID 

constant  for  t  >  0  at  one  surface  (x  »  0)  and  insulated  at  the  other 
(x  *  i).  This  represents  the  worst  heating  condition  possible  since 
the  heat  at  x  ■  i  will  actually  be  transferred  away  by  radiation  and 
convection  to  the  surrounds.  The  conduction  of  heat  through  the  slab 
is  governed  by  the  well  known  equation 


k  o  T  dT 

wc  >2  57 


The  boundary  conditions  are 


T(o,t) 


I  ('.=)  .  0 


Equation  (1)  can  be  easily  solved  by  Laplace  Transform  techniques.  Hence, 
defining  D  ■  k/wc  -  thermal  diffusivity,  Eq.  (1)  in  the  transform  T(x,s) 


becomes 


d2T  - 

D  =-4  «  s  T 

dx 


which  has  the  solution 


_  -  V  s/D  x  V  s/D  x 

T(x,s)  -  A  e  +  B  e 


The  transformed  boundary  conditions  are 


T(o,s)  -  f 


106 


£2  (i,s)  _  Q 

dx 


(6) 


Substitution  of  (5)  and  (6)  into  (4)  yields  two  simultaneous  equations 
which  can  be  readily  solved  for  the  arbitrary  constants  A  and  B. 
With  A  and  B  known  the  solution  is 


T  (x.s)  ^  cosh ■>/a/D~V  (l-x) 


m 


8  cosh  •J s/d  i 


(7) 


For  inversion  purposes,  express  the  hyperbolic  functions  in  the 
exponential  form,  and  expand  in  a  series  by  the  binomial  theorem. 
Hence  (7)  gives 


-  Vs/D  x  -  Vs/D  (21-x) 
1  e  +  e 


T  s 
m 


1  +  e-2^7r, 

.  [“  -  V  s/d'  x  -  'Js/D  (21-x)  1 

s  Le  +  e  J 


-  </  s/D  x  -  'J  s/D  (21-x)  ~|  »  -2n  '/s/d"  1 

e  +  e  I  Z  (-l)n  e 

n-0 


Hence 


-  ,  oo  -  >/s/D(2nl+x)  oo 

=  -  2  (-l)n  e  +  -  2  (-1)“  e 

T  s  .  '  s  -  '  ' 

m  n=0  n*0 


This  can  readily  be  inverted  since  from  ref.  32 


[  (2n+2)l  -  x] 

(8) 


-IT.  Ws/DxT 

]j  •  J  ■ 


erfc 


Jot 


(9) 


The  notation  (erfc)  is  given  by 


erfc  x  ■  1  -  erf  x 


where  (erf)  is  the  error  function.  The  inversion  of  (8)  is  readily 
obtained  by  the  use  of  (9),  and  hence  the  solution  for  the  temperature 
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In  Che  coaclng  is 

-  £  (-l)n  erfc  2b£JLx  +  £  (-l)n  erfc  -  -  (10) 

m  n-0  2  vDt  n-0  2  VDt 

This  solution  converges  quickly,  and  in  order  to  Illustrate  the  time 
scales  for  the  coating  to  experience  a  significant  temperature  rise  it 
will  be  necessary  to  evaluate  (10)  for  some  specific  times  and  distances, 
ttie  properties  of  the  coating  are  as  follows. 


k  »  0,00057  cal/cm  sec  °C 
7  -  1.239  gms/cc 

c  -  0.27  cal/gm  °C 


k_ 

7c 


0.00057 


(1.289)(0.27) 

0.030  inches  -  0.076  cm. 


0.0017  cm  /sec. 


*  0.294/sec. 

£ 

At  x  *  0,  obviously  T  =  Tm>  since  the  coating  is  separated  by  a  glue 
interface  of  the  same  properties  as  the  coating.  We  will  assume  the  glue 
to  be  0.005  inches  thick  and  calculate  the  temperature  history  of  x/J  - 
0.167  from 


m 


£ 

n-0 


(-l)n  erfc 


n  +  0,084 
0.54  */t 


+  £ 
n-0 


(-l)n  erfc 


n  +  0.92 
0.54  >ft 


This  calculation  results  in  the  following  temperature  history  at  x  -  0.005 


t 


inches 
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t(secs.)  0  0.010  0.050  .  0.100  0.500  1.000 

0  0.038  0.32  0.48  0.77  0.84 

T 

m 

As  a  result  of  this  calculation,  it  can  be  concluded  that  the  coating 

remains  at  substantially  its  initial  temperature  during  the  first  ten 

milliseconds  after  bank  discharge.  The  significance  of  this  result  lies 

in  the  fact  that  the  optical  properties  of  the  coating  remain  constant, 

with  the  exception  of  possible  dynamic  effects,  since  all  observations 

were  made  within  about  the  first  300  >isecs.  In  order  to  verify  the 

validity  of  this  solution,  a  literature  search  revealed  that  this  result 

* 

is  identical  with  that  previously  obtained  by  Carslaw  and  Jaeger  for  a 
very  similar  problem.  These  results  are  conservative  since  it  was  assumed 
that  the  coating  was  insulated  over  its  outer  surface  and  that  the 
temperature  of  the  bar  remained  constant.  This  latter  assumption  has 
been  investigated  in  the  previous  section,  where  it  was  shown  that  the 
temperature  remains  nearly  constant  with  respect  to  a  millisecond  time 
scale. 


4 


If 

H.  S.  Carslaw  and  J.  C.  Jaeger, 
Press,  1959,  2nd  Ed.,  pp  96-99. 


Conduction  of  Heat  in  Solids. 
(Note  Eq.  9,  p.  97). 


Oxford 
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4.  Evaluation  of  the  Secondary  Principal  Stress 

In  order  to  simplify  the  calculation  of  the  axial  stress  from  the 
stress'optic  law  (3.11),  the  secondary  principal  stress  was  assumed  to 
be  negligible.  This  simplifying  assumption  is  justified  as  follows. 
The  secondary  principal  stress  0^  arises  from  the  inertial  loading 
of  the  coating  which  must  be  accelerated  in  the  y-direction  due  to  the 
lateral  expansion  and  contraction  of  the  bar.  Hence,  the  inertial 
stress  0^  in  the  coating  is  (refer  to  Fig.  1) 


where  m^  -  mass  of  coating 
•  • 

y  ■  lateral  acceleration 
h  «*  coating  thickness 
L  ■  length  of  coating 
The  lateral  acceleration  is  given  by 

«  t  • 

y  **  Va  e 

x 

where  a  =  the  sonic  velocity,  and  6^  ■  axial  strain  rate. 

The  axial  strain  rate  can  be  determined  from  the  theoretical  solutions 
of  Sec.  2.  Hence,  assuming  Case  II  heating,  Eq.  (2.34)  gives 

e  -  30  in/ in/sec. 
max 

and  since 

m  <■  p  n  6  L 
c  c 


the  maximum  stress  is 
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6  v-  *  '* 


Substituting  some  representative  magnitudes. 


„  (0,045) (0.050) (0,30) (0.2  x  10°)(30)  .  „  . 

y  386  U  PS1 

1  max 


Hence,  o  is  very  small  compared  to  the  axial  stress  o  ,  and  can, 

y  x 

therefore,  be  neglected. 


i 


t 


I 


4 
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5.  Reinforcing  Effect  of  the  Photoelastic  Coating 

In  order  to  determine  the  reinforcing  effect  of  the  coating  consider 
a  rectangular  bar  of  uniform  section  A  and  density  to  be  coated  uni- 
formly  on  two  of  its  sides  with  a  material  having  mass  density  Pc  and 
cross-sectional  area  A£  as  shown  in  Fig.  6A.  With  the  passage  of  a 
longitudinal  wave  along  the  bar,  the  elemental  length  PQ  will  be 
displaced  to  the  new  position  P’Q*.  If  the  coating  thickness  is  small 
relative  to  the  bar  thickness,  then  we  can  assume  that  shear  strains  in 
the  coating  are  negligible  relative  to  the  axial  strain.  Further,  we 
will  assume  that  the  axial  strain  in  the  coating  is  equal  to  the  uniaxial 
strain  in  the  bar,  i.e. 


e 


c 


e 

s 


auc  aus  du 

sr  -  sr 


(5.1) 


The  influence  of  the  coating  on  the  motion  of  the  bar  can  be  determined 
by  considering  the  forces  acting  on  the  element  P'Q'.  A  Force  balance 
on  these  elements  yields. 

da  d2u 

Bar:  -  2F  +  Afi  ^  dx  -  A^  dx  (5.2) 


dcc 

Coating:  F  -  Ac  •  dx 


A  p 
c  c 


(5.3) 


Eliminating  F  between  Eqs.  (5.2)  and  (5.3),  and  substituting  (5.1)  we 

get, 


sr  Afl  sr  r  +  A8p8/;p 


dfu 

dt2 


(5.4) 


Differentiation  of  (5.4)  with  respect  to  x  yields 
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b2a  A  S2o 

_ 8  C  _ C 


(5.5) 


The  stress-strain  relations  are 


r  +  as T(t) 

s 


(5.6) 


(5.7) 


3x2 


(5.8) 


Hence,  substitution  of  (5.6)  and  (5.8)  into  (5.5)  yields  the  equation 


d2ff  tn  d2g  d2T 

— r-  m  — -  +  m  p  a  — r 

a*2  a2  at2  8  8  at2 


(5.9) 


where 


A  p 

1+JL  JL 

As  Ps 
'As  Es 


(5.10) 


If  we  define 


a'  - 


-  £ 


(5.11) 


"/in'  t 


(5.12) 


Hence,  the  corrected  equation  (2.8)  is 


i  a fa 

a ,2  at2 


+  m  p  a 
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which  has  the  same  solution  given  in  Sec.  2  except  that  t  is  exchanged 
for  t'.  It  should  be  noted  that  m  appears  in  the  terms  multiplying 
the  series  solution  but  cancels  out.  For  example,  in  Eq.  (2.26) 

2 

paa2T  «  m  p  a  a 1 2  T  »  mpa—  =  ESI  =  a 
m  m  m  m  bt 

Hence,  only  the  period  x  is  affected.  Consequently  the  stress  amplitude 
is  increased  slightly  due  to  the  dependence  upon  the  ratio  t./r  which 
decreases  as  a  result  of  the  increase  of  x  to  t'  by  the  factor  •Jn 7. 

This  correction  is  not  exact  since  the  strain  at  the  free  surface  of 

the  coating  is  different  from  that  at  the  interface  and  shear-strains  are 

undoubtedly  present.  An  exact  solution  of  this  problem  would  be  extremely 

difficult  to  obtain,  particularly  in  view  of  the  fact  that  inertia  effects 

must  be  included.  The  exact  solution  for  static  loading  has  been  worked 
40 

out  by  Duffy  who  has  shown  that  the  birefringence  observed  perpendicular 
to  the  surface  is  significantly  dependent  upon  the  strain  gradient  and 
curvature  at  the  metal  surface.  Since  the  birefringence  is  observed 
parallel  to  the  surface  in  this  work,  these  results  do  not  directly 
apply  here.  This  analysis  illustrates  nicely  the  nature  of  the  required 
exact  solution  which  would  be  further  complicated  by  the  addition  of 
inertial  terms. 

While  the  correction  factor  (5,10)  is  not  exact,  it  is  considered 
sufficiently  accurate  for  engineering  purposes.  Since  the  derivation  is 
based  on  the  assumption  that  the  bar  is  coated  over  its  entire  length 
rather  than  only  over  a  portion,  this  correction  is  undoubtedly  conserva¬ 
tive.  This  is  apparent  from  an  energy  consideration  since  in  most  cases 
the  total  mass  of  the  specimen  is  much  lower  than  that  of  the  coating. 

As  pointed  out  in  Sec.  5,  two  of  the  specimens  were  tested  with  different 
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coating  lengths.  It  was  observed  that  some  coating  effect  on  the  response 
may  have  been  present,  but  only  when  the  coating  and  specimen  volumes  were 
nearly  equal.  Further,  the  period  of  oscillation  measured  experimentally 
is  larger  than  the  theoretical  period  t  by  an  amount  very  closely 
approximated  by  (5.12).  Also,  the  physical  nature  of  (5.10)  is  reasonable 
since  it  indicates  that  the  inertia  effect  appears  in  the  numerator  while 
the  stiffness  appears  in  the  denominator  which  is  consistent  with  the 
definition  of  natural  period  of  oscillation.  Hence,  it  is  evident  that 
confidence  in  the  validity  of  the  correction  factor  as  derived  is 


justifiable. 
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6.  Calibration  of  the  Photoelastic  Coating 

In  order  to  determine  if  the  strain-optic  coefficient  of  the  coating 
is  rate  dependent,  it  was  necessary  to  perform  a  dynamic  calibration. 

Since  there  is  no  data  available  in  the  literature  on  the  specific  bire- 
fringent  material  used  in  this  work,  the  calibration  must  be  carried  out 
at  substantially  the  same  loading  rates  experienced  in  the  experiment. 

Also,  the  dynamic  strain  must  be  measured  independently  as  well  as 
simultaneously  with  the  photoelastic  method.  Since  the  effects  of 
rapidly  heating  an  elastic  rod  are  identical  to  the  effects  of  heating 
a  rod  slowly,  compressing  it  back  to  its  original  length,  and  then 
releasing  it  suddenly,  it  was  apparent  that  an  effective  calibration 
could  be  performed  in  the  following  manner. 

Rather  than  precompress  a  rod  and  design  a  complicated  device  to 
provide  a  sudden  release,  it  was  decided  to  prestretch  a  vertical  rod 
clamped  at  its  upper  end,  by  loading  it  with  a  weight  suspended  from 
a  thin  wire  fastened  at  the  free  end  of  the  rod.  Sudden  release  of  the 
weight  was  accomplished  by  utilizing  the  capacitor  band  to  explote  the 
lower  portion  of  the  wire.  Hence,  the  0.125  x  0.250  x  5.00  inch 
specimen  (run  2)  was  clamped  in  the  heating  jig  as  shown  in  Fig.  7A  and 
preloaded  with  an  80  lb.  weight  suspended  by  a  0.016  inch  diameter  music 
wire  approximately  15  inches  long.  Two  350  ohm  semi-conductor  strain 
gages  were  mounted  on  the  uncoated  sides  of  the  bar  at  the  same  location 
as  the  coating.  These  gages  were  connected  in  series  to  cancel  any 
effects  due  to  lateral  bending  and  insure  measurement  of  only  longitudinal 
strain.  With  the  gages  in  a  simple  voltage  divider  circuit  as  shown,  their 
output  could  be  read  directly  on  the  535  Tektronix  oscilloscope.  Since 
these  gages  have  a  gage  factor  of  123,  small  strains  are  easily  measured 


fig.  7a.  Dynamic  calibration  of  photoelastic  coating. 
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with  little  amplification.  The  strain  was  also  measured  independently  with 
the  photoelastic  coating  technique  described  in  Sec.  3  by  using  all  existing 
photographic  equipment. 

The  sequence  of  events  was  very  similar  to  the  rapid  heating  tests. 

The  room  was  darkened,  camera  brought  up  to  speed,  and  with  closure  of  the 
"fire"  switch  the  synchronizing  pulse  triggered  the  oscilloscopes,  dis¬ 
charged  the  capacitor  bank  to  explode  the  wire,  and  after  a  suitable  delay 
discharged  the  light  source  and  bridge-wire.  In  this  manner,  the  subsequent 
longitudinal  oscillations  of  the  rod  were  observed  independently  and 
simultaneously  by  electronic  and  photoelastic  methods.  The  advantages  of 
this  method  are  that  it  was  simple  and  provided  a  test  at  a  representative 
frequency  -  in  this  case,  100  nsec.  period  since  the  rod  was  5  inches  long. 
Also,  sudden  release  was  obtained  as  the  mass  of  the  wire  is  small  relative 
to  the  bar  mass.  Since  the  unloading  wave  in  the  wire  travels  at  the  sonic 
velocity,  the  length  of  the  wire  was  selected  so  that  the  transit  time 
(L/a  ■  75  jxsec.)  was  greater  than  the  discharge  time  of  the  capacitor  bank. 
Hence,  the  strain  gage  output  was  obtained  free  of  electromagnetic  noise 
from  the  exploding  wire.  The  light  source  delay  was  therefore  selected  to 
be  approximately  equal  to  this  transit  time. 

The  M601  film  strip  provided  the  fringe  movements  from  which  the  fringe 
order  was  obtained  at  a  point.  In  order  to  calculate  strain,  the  static 
value  of  the  strain-optic  coefficient  was  obtained  by  loading  the  bar  slowly 
In  a  standard  tensile  testing  machine  and  observing  the  fringe  movements 
with  a  telescope.  Also,  the  output  of  the  strain  gages  was  observed.  It 


was  found  that 
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and 

e2  -  1.09  eQ  (6.2) 

where  eQ  is  the  output  of  the  voltage  divider  in  millivolts,  €  is  the 
strain  in  microinches/in.,  and  h  is  the  width  of  the  coating  (0.250 
inches).  Using  (6.1),  the  dynamic  strain  was  determined  from  the  film 
record  in  the  same  manner  as  previously  described.  This  result  is 
compared  with  the  strain  gage  record  in  Fig.  8A. 

Since  the  static  strain-optic  coefficient  was  used,  it  is  evident 
from  Fig.  8A  that  good  agreement  exists  and  the  dynamic  strain-optic 
coefficient  is  not  seriously  reduced  at  this  rate  of  loading.  Hence, 
it  can  be  concluded  that  the  ratio  of  dynamic  to  static  strain-optic 
coefficients  can  be  taken  as  unity  for  thi6  work  and  the  simplified 
equation  (3.18)  used  for  data  reduction. 

This  calibration  technique  is  extremely  useful  since  it  can  provide 
tests  over  a  range  of  frequencies  restricted  only  by  the -length  of  the 
rod.  The  only  difficulty  arises  in  selecting  the  weight.  In  order  to 
provide  enough  space  for  the  strain  gage  mounting,  the  cross-sectional 
area  of  the  bar  cannot  be  much  less  than  the  one  used.  Since  the  wire 
diameter  should  be  small,  the  weight  magnitude  is  limited,  which  in  turn 
limits  the  initial  axial  strain  present  in  the  bar  during  the  test. 
Consequently,  the  small  fringe  movements  observed  in  the  coating  are 
sometimes  difficult  to  measure  visually.  The  data  shown  in  Fig.  8A 
was  obtained  from  measurement  of  fringe  orders  which  did  not  exceed  0.14. 
Hence,  there  was  a  fair  amount  of  scatter  present,  and  it  was  necessary 
to  pick  those  points  which  fell  on  the  "best  curve"  drawn  through  all  of 
the  experimental  points.  However,  this  is  not  a  serious  compromise  since 
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care  was  taken  to  select  only  those  points  which  represented  the  average 
trend  of  the  data  so  that  there  were  as  many  points  above,  as  well  as  below 
the  strain  gage  record.  The  fact  that  the  period  of  oscillation  was  also 
reproduced  accurately  provides  sufficient  evidence  to  conclude  that  dynamic 
optical  effects  can  be  neglected. 
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FIG.  8A.  COMPARISON  OF  PHOTOELASTIC  COATING  DATA  WITH  STRAIN  GAGE  DATA 


